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Abstract 

This  thesis  proposes  a  new  direct-sequence  spread  spectrum  (DS/SS)  code  phase  tracking 
loop  which  mitigates  the  effects  of  multipath  interference  on  code  phase  tracking  error;  such  errors 
can  translate  to  significant  range  measurement  errors  in  DS/SS  ranging  systems  such  as  Global 
Positioning  System  (GPS).  The  new  code  tracking  loop,  called  the  modified  RAKE  delay-lock  loop 
(MRDLL),  uses  mcodmum-likelihood  (ML)  signal  parameter  estimation  to  determine  the  amplitude, 
carrier  phase,  and  relative  propagation  delay  of  both  a  direct-path  and  a  reflected  signal;  a  multiple- 
correlator  code  phase  tracking  loop  then  exploits  these  ML  signal  estimates  to  remove  the  tracking 
error  introduced  by  the  reflection.  A  preliminary  analysis  showed  that  the  MRDLL’s  linear  tracking 
region  varied  with  the  reflected  signal  parameters;  therefore,  an  adaptive  loop  controller  (ALC)  was 
introduced  to  allow  the  loop  designer  to  fix  dynamic  specifications  such  as  loop  natural  frequency. 
Analysis  and  computer  simulations  demonstrated  that,  when  multipath  was  present,  the  MRDLL 
exhibited  a  significantly  lower  steady-state  code  phase  tracking  error  than  that  of  the  standard 
non-coherent  delay-lock  loop  (NCDLL),  which  is  typically  used  in  GPS  receivers.  In  an  ideal 
multipath-free  environment,  the  NCDLL  is  still  the  best  choice  for  code  phase  tracking. 

This  GPS  receiver  design  technology  will  benefit  the  entire  aviation  community  by  eliminating 
or  reducing  the  dominant  source  of  error  in  differential  GPS-based  instrument  landing  systems, 
resulting  in  improved  ILS  safety,  reliability,  and  integrity.  In  relative  GPS  applications,  such 
as  precision  guided  munitions,  lethality  is  improved  via  elimination  of  multipath  contribution  to 
targeting  error.  In  all  GPS  applications,  the  (no  longer  dominant)  error  contribution  of  multipath 
can  be  eliminated,  yielding  enhanced  positioning  accuracy. 
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ANALYSIS  AND  SIMULATION  OF  A  NEW  CODE  TRACKING  LOOP 
FOR  GPS  MULTIPATH  MITIGATION 


I.  Introduction 

1.1  Overview 

The  satellite-based  NAVSTAR  Global  Positioning  System  (GPS)  has  ushered  in  a  new  era  in 
world-wide  precision  navigation.  Both  the  military  and  civilian  communities  use  GPS  in  a  variety 
of  applications  as  users  rely  on  GPS  satellite  signals  to  estimate  current  latitude,  longitude,  and 
elevation.  The  United  States  Air  Force  (USAF),  for  example,  relies  on  GPS  receivers  to  provide 
accurate  position  data  for  many  of  its  aircraft  navigation  systems. 

In  the  ongoing  search  for  greater  navigation  accuracy,  research  efforts  continue  to  focus  on 
understanding  and  reducing  the  effects  of  real-time  errors  at  the  GPS  receiver.  Three  significant 
error  sources  are  satellite  and  receiver  clock  biases,  atmospheric  effects  on  the  GPS  signal,  and 
multipath  signal  reflections.  Two  of  these  error  sources  (clock  bias  and  atmospheric  effect)  can 
often  be  significantly  reduced  by  using  differencing  techniques  such  as  differential  GPS  (DGPS). 
Multipath  effects,  however,  generally  cannot  be  reduced  by  differencing  techniques;  therefore,  re¬ 
ducing  the  effects  of  multipath  interference  on  the  GPS  receiver  has  received  much  attention.  This 
thesis  examines  a  proposed  GPS  receiver  code  tracking  loop  designed  to  mitigate  the  effects  of 
multipath  interference. 

1.2  Background  Information 

1.2.1  The  Multipath  Signal.  Multipath  signals  occur  when  the  original  GPS  line-of-sight 
(LOS)  satellite  signal  reflects  off  surrounding  objects  or  the  ground  plane  (Figure  1).  At  the  GPS 
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Figure  1.  Multipath  signal  reflection, 

receiver,  this  reflection  produces  an  attenuated  replica  of  the  LOS  signal  which  is  delayed  due  to 
the  longer  propagation  paths. 

The  general  multipath  signal  can  be  expressed  as  the  sum  of  several  time-shifted  and  atten¬ 
uated  versions  of  the  original  transmitted  signal: 


^  ^  O-jS  (i  Ti) 


where  s(t)  is  the  original  transmitted  signal,  L  is  the  number  of  reflected  paths  present,  is  the 
attenuation  coefficient  of  the  ^th  signal,  and  is  the  propagation  delay  of  the  zth  signal  (i  =  0 
corresponds  to  the  direct  path  LOS  signal). 

For  a  single  GPS  satellite,  the  transmitted  signal  can  be  represented  as 


s(^t)  =  c{t)m{t)\/^  cos{u)ct  +  <t>) 


where  c{t)  is  a  pseudo-random  noise  (PN)  direct-sequence  spread-spectrum  (DS/SS)  code,  m{t) 
represents  binary  phase-shift  keyed  (BPSK)  data  modulation  at  50  bits  per  second  (referred  to 
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propagation  delay 

Figure  2,  Average  power-delay  profile  for  a  multipath  signal  with  a  single  reflection. 

as  the  ‘navigation  message’),  P  is  the  signal  power  in  Watts,  and  (f>  is  the  initial  transmitted 
phase  of  the  signal  in  radians.  GPS  transmits  at  two  different  carrier  frequencies,  Wc=  27r fc  (in 
radians/sec);  the  first  carrier,  designated  LI,  is  at  fc  =  1575.42  MHz  while  the  second  carrier,  L2, 
is  at  fc  =  1227.6  MHz  (4). 

When  modeling  the  attenuation  coefficients,  a^,  it  is  important  to  consider  both  their  dis¬ 
tribution  and  their  values  relative  to  each  other.  There  are  two  scenarios  to  be  considered  when 
modeling  the  amplitude  distribution  of  the  received  signal  (11): 

1.  A  single  reflection  is  present.  In  this  case,  L  =  1  in  Equation  1.  A  typical  power-delay 
profile  for  ao  and  ai  is  shown  in  Figure  2.  This  scenario  is  valid  when  there  is  one  dominant 
reflector  present  such  as  a  large  building  or  body  of  water.  Often,  ao  and  ai  are  assumed  to 
be  uniformly  distributed  between  0  and  1. 

2.  Multiple  reflections  are  present.  In  this  case,  L  >  1  in  Equation  1,  and  the  a^’s  are 
distributed  according  to  a  Rayleigh  distribution.  This  scenario  is  valid  when  many  reflectors 
are  present,  such  as  buildings  or  trees.  A  typical  power-delay  profile  for  the  multiple-reflection 
case  is  shown  in  Figure  3 
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Figure  3.  Average  power-delay  profile  for  a  multipath  signal  with  multiple  reflections. 

Notice  that  the  amplitudes  of  the  reflected  signals  are  generally  less  than  the  direct  path 
signal.  According  to  (11),  this  condition  will  occur  for  the  following  reasons: 


1.  The  reflected  signal  travels  a  greater  distance.  Therefore,  it  experiences  a  greater  free 
space  loss.  However,  this  attenuation  has  a  minimal  effect  on  GPS  receiver  ranging  errors 
because,  as  will  be  seen,  only  reflections  with  a  propagation  delay  just  slightly  greater  than 
that  of  the  direct-path  signal  have  a  significant  impact. 

2.  The  reflector  surface  causes  attenuation.  The  amount  of  this  attenuation  is  dependent 
upon  the  material  of  the  reflector  and  the  incident  angle  of  the  reflection.  For  very  low 
incident  angles,  the  attenuation  is  negligible. 

3.  The  receiver  antenna  causes  attenuation.  This  is  due  to  the  antenna’s  gain  pattern 
and  the  orthogonal  polarization  between  the  received  signals.  Generally,  in  the  case  of  GPS, 
a  single  reflection  has  a  left-hand  circular  polarization  while  the  direct  path  signal  has  a 
right-hand  polarization.  Therefore,  the  two  received  signal  components  are  orthogonal  to  one 
another. 
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1.2.2  The  GPS  Spreading  Code.  In  order  to  understand  the  operation  of  the  GPS  receiver 
and  the  effects  of  multipath,  it  is  first  necessary  to  examine  the  properties  of  the  DS/SS  code,  c{t). 
This  code  is  a  BPSK  code  generated  at  a  very  high  frequency  relative  to  the  data  message  rate. 
To  distinguish  the  spreading  code  from  the  data  message,  each  code  pulse  of  is  called  a  ‘chip’  (as 
opposed  to  a  ‘bit’)  and  the  code  is  said  to  have  a  ‘chip  period’  of  Tc  seconds.  Furthermore,  the 
spreading  code  itself  is  periodic  with  period  NTc  seconds,  where  N  is  the  number  of  chips  in  a 
period.  The  actual  values  of  Tc  and  N  differ  depending  on  whether  the  user  is  receiving  the  lower- 
rate  coarse/acquisition  (C/A)  code  or  the  higher-rate  precision  code  (P-code).  The  properties  of 
these  two  codes  are  summarized  in  Table  1  which  is  based  on  information  presented  in  (4). 


Table  1.  GPS  spreading  code  parameters. 


Parameter 

C/A  code 

P-code 

chip  rate  =  1/Tc 

1.023  MHz 

10.23  MHz 

code  period  =  N 

1023  chips;  (1ms) 

«  6  X  10^^  chips;  (Iweek) 

carrier  band 

Ll 

Ll,  L2 

At  the  transmitter,  mixing  c{t)  with  the  data-modulated  carrier  has  the  effect  of  ‘spreading’ 
the  GPS  signal’s  bandwidth  and  significantly  reducing  the  effective  received  signal-to-noise  ratio 
(SNR).  A  receiver  can  recover  the  transmitted  signal  by  mixing  the  received  signal  with  an  internally 
generated  and  synchronized  replica  code;  this  is  known  as  ‘despreading’  the  signal.  The  replica  code 
is  generated  by  a  code  tracking  loop  in  the  receiver  which  is  responsible  for  tracking  changes  in  the 
propagation  delay  (or  code  phase)  of  the  incoming  DS/SS  code. 

An  important  property  of  any  DS/SS  code  is  its  autocorrelation  function,  Rc{^)t  which  is 
defined  as 

Rc(n)  =  —  jf  c{t)c{t  -t-  Q.Tc)dt  .  (3) 
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Figure  4.  DS/SS  spreading  code  autocorrelation  function. 

This  function  is  illustrated  in  Figure  4  for  a  maximal-length  PN  code  (similar  to  the  Gold  codes 
used  by  GPS);  it  is  periodic  with  period  N  and  can  be  written  as  (5) 

l-Q(l  +  i)  0<fi<l 

i2c(0)=<  l<n<{N-l)  ■  (4) 

[n - (iv - 1)] (1  +  i)  - i  {N-i)<n<N 

As  can  be  seen  in  Table  1,  iV  >  1  for  the  GPS  signal.  Therefore,  the  code  autocorrelation 
function  can  be  approximated  as 


1  -  N  |fi|  <  1 

0  elsewhere 


(5) 


L£,S  The  code  tracking  loop.  As  mentioned  before  in  Section  1.2.2,  the  GPS  receiver’s 
code  tracking  loop  tracks  the  code  phase  of  the  incoming  DS/SS  signal,  c(t).  Tracking  the  code 
phase  serves  two  main  purposes: 

1.  Tracking  provides  an  estimate  of  the  magnitude  of  time  shift  required  to  maximize  the  cor¬ 
relation  between  the  incoming  signal  and  the  receiver’s  internally  generated  ‘on-time’  code; 
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this  information  is  used  by  the  receiver  to  calculate  an  initial  user-to-satellite  range  estimate, 
known  as  the  pseudorange  (PR)  measurement. 

2.  The  synchronized  replica  signal  is  used  to  despread  the  GPS  signal;  this  despread  signal  is 
then  passed  to  the  receiver’s  carrier  tracking  loop  for  demodulation  of  the  data  message, 
m{t).  The  data  message  contains  satellite  information  that  is  needed  to  enable  the  receiver 
to  compute  its  position,  velocity,  and  clock  bias. 

The  typical  GPS  receiver  uses  the  standard  non-coherent  delay-lock-loop  (NCDLL)  as  its  code 
tracking  loop.  This  loop  is  often  called  the  ‘early-late’  delay-lock  loop,  referring  to  the  advanced  and 
delayed  versions  of  the  code  replica  generated  by  the  voltage-controlled  clock  (VCC).  Furthermore, 
the  relative  delay  spacing,  ATc,  between  the  two  generated  codes  is  called  the  ‘early-late  spacing’ 
(typically,  A  =  1.0),  The  the  theory  of  operation  and  an  analysis  of  the  NCDLL  will  be  presented 
in  Chapter  11. 

1.2,4  Effects  of  multipath  on  the  code  tracking  loop.  The  code  tracking  loop’s  mission  is 
to  accurately  track  the  code  phase  of  the  direct  path  LOS  signal.  However,  when  multipath  signals 
arrive  at  the  GPS  receiver,  they  introduce  a  phase  tracking  error  in  the  NCDLL.  The  significance 
of  these  tracking  errors  is  larger  than  one  might  initially  think.  For  instance,  for  C/A  code,  a  code 
phase  tracking  error,  e,  of  one-tenth  of  a  chip  period,  T^,  gives: 

€  =  0.17;«98ns,  (6) 

which  may  seem  rather  small.  However,  representing  this  tracking  error  in  terms  of  distance  (e^) 
by  multiplying  by  the  speed  of  light,  c,  gives 

ed^ce^  {3  X  10®  m/s)  (9.8  x  10®  s)  =  29,4  m.  (7) 
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Therefore,  a  code  tracking  error  of  only  98  ns  corresponds  to  a  range  measurement  error  of  almost 
30  m!  It  has  been  shown  (11)  that  multipath  tracking  errors  are  significant  only  when  the  relative 
delay  between  the  multipath  and  direct-path  signals  is  less  than  1.5Tc;  this  is  due  to  the  inherent 
multipath  resistance  of  DS/SS  spreading  codes  and  is  based  on  properties  of  the  code  autocorrela¬ 
tion  function  defined  in  Equation  3,  The  effects  of  multipath  on  the  NCDLL’s  code  phase  tracking 
performance  is  examined  in  detail  in  Chapter  11. 

1.2,5  Proposed  Mitigation  Techniques.  A  significant  amount  of  research  has  gone  into  the 
mitigation  of  GPS  multipath  effects.  This  section  looks  at  three  mitigation  techniques  that  have 
been  proposed  by  various  authors: 

1.  Narrow  correlator  spacing 

2.  Minimum  mean-square  error  (MMSE)  estimation 

3.  Multiple-correlator  DLL 

1.2.5. 1  Narrow  Correlator  Spacing.  It  was  shown  by  Van  Dierendonck,  Fenton,  and 
Ford  that  reducing  the  early-late  spacing  in  the  NCDLL  could  reduce  GPS  multipath  errors  under 
certain  conditions  (10).  Experiments  performed  with  a  C/A  code  GPS  receiver  demonstrated  that 
reducing  the  standard  1.0-chip  spacing  to  0.1-chip  spacing  reduced  the  tracking  error  caused  by 
multipath  interference.  However,  the  authors  in  (10)  admit  that  the  problem  with  the  narrow 
correlator  spacing  technique  is  that  it  currently  can  only  be  used  with  C/A  code,  which  is  intended 
primarily  for  civilian  use.  The  NCDLL  with  narrow  spacing  cannot  accurately  track  the  higher- 
frequency  P-code,  which  is  used  in  important  military  applications. 

1.2. 5. 2  MMSE  Estimation.  A  second  multipath  mitigation  approach  proposed 
replacing  the  NCDLL  altogether  with  an  alternative  direct-path  delay  estimator.  This  alternative 
estimator  was  derived  by  Weill  as  the  minimum  mean-square  error  (MMSE)  estimator  for  direct 
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path  delay  (14).  It  had  already  been  shown  that  the  NCDLL  approached  the  performance  of  the 
maximum-likelihood  estimator  (MLE)  for  one-path  signals  (13).  However,  while  the  MLE  was 
optimum  for  the  one-path  case,  Weill  showed  in  (14)  that  the  MMSE  estimator  performed  better 
in  the  presence  of  multipath  because  no  other  estimator  exhibited  a  lower  root  mean-square  (rms) 
estimation  error;  the  only  exception  to  this  occurred  when  multipath/LOS  path  separations  were 
extremely  small  (approximately  1-2  meters).  The  main  disadvantage  of  the  MMSE  estimator  lies 
in  implementation.  It  is  currently  difficult  to  implement  the  true  MMSE  estimator  for  direct-path 
delay,  because  the  MMSE  requires  two  quadruple  integrations  which  are  very  computationally 
intensive. 


1.2. 5. 3  Multiple-correlator  estimation.  A  third  multipath  mitigation  scheme  is  the 
multiple  correlator  estimation  (MCE).  As  illustrated  in  Figure  5,  the  received  signal  is  correlated 
with  many  delayed  code  replicas  (c(f  -  f  -  PkTc)  in  Figure  5)  provided  by  the  code  tracking  loop. 
As  will  be  seen  in  Chapter  III,  this  effectively  ‘samples’  the  code  autocorrelation  function,  Rc{^), 
at  fl  =  the  resulting  correlator  outputs,  Rk,  are  sampled  and  digitally  processed  to  estimate  the 
gain  and  phase  of  the  direct-path  and  reflected  signal  components.  These  estimates  are  then  used 
by  the  code  tracking  loop  to  reduce  any  code  tracking  errors  introduced  by  multipath  interference. 
One  example  of  the  MCE  approach  is  the  Multipath  Estimating  Delay  Lock  Loop  (MEDLL).  The 
MEDLL  uses  MCE  to  remove  the  estimated  multipath  signal  components  and  feed  the  estimated 
direct-path  signal  to  a  standard  NCDLL  (9). 

The  MCE  approach  was  also  used  by  Sheen  and  Stuber  in  their  proposed  RAKE  delay- 
lock  loop  (RDLL)  (7)  (8).  This  approach  differed  from  the  MEDLL  in  that  the  NCDLL  was  not 
used;  instead,  the  code  tracking  loop  was  completely  redesigned.  The  RDLL  featured  a  relatively 
simple  multiple-correlator  loop  design  that  demonstrated  a  significant  improvement  in  tracking  error 
performance  over  the  standard  NCDLL  in  the  presence  of  multi-user  interference.  However,  unlike 
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Figure  5.  Multiple-correlator  estimation. 

the  MEDLL,  the  RDLL  Wcis  not  designed  for  GPS  applications,  but  for  code-division-multiple-access 
(CDMA)  mobile  radio  communications. 

1.3  Problem  Statement 

The  RDLL  was  designed  to  operate  in  a  multiple-user  CDMA  environment  with  relative 
delay  spacings  that  were  integer  multiples  of  the  chip  period,  Tc-  As  stated  in  Section  1.2,4,  GPS 
code  tracking  loop  errors  due  to  multipath  occur  only  for  multipath  delay  spacings  that  are  very 
small  (i.e.,  fractional  multiples  of  Tc).  Therefore,  for  the  RDLL  to  be  effective  in  GPS  multipath 
environments,  it  must  be  modified  to  mitigate  the  effects  of  multipath  signals  with  small  delays 
relative  to  the  direct-path  signal. 

1.4  Thesis  Objectives 

This  thesis  proposes  a  GPS  receiver  code  tracking  loop,  the  modified  RAKE  delay-lock  loop 
(MRDLL),  designed  to  mitigate  the  effects  of  multipath  interference.  This  new  tracking  loop  is 
a  modified  version  of  the  RDLL  introduced  by  Sheen  and  Stuber.  The  proposed  loop  includes  a 
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multiple-correlator  estimation  unit  (MCEU)  and  introduces  an  original  adaptive  loop  controller 
(ALC)  not  present  in  the  RDLL. 

The  main  objectives  of  this  research  were  to: 

1.  Verify  the  effects  of  multipath  on  NCDLL  code  tracking  performance. 

2.  Adapt  the  existing  RDLL  design  for  solution  of  the  GPS  multipath  problem. 

3.  Design  an  original  multipath  gain/phase  estimation  scheme  for  use  with  the  MRDLL. 

4.  Compare  predicted  and  simulated  performance  of  the  MRDLL  to  that  of  the  NCDLL  in  the 
presence  of  multipath. 

5.  Investigate  performance  of  the  MRDLL  in  the  presence  of  additive  white  Gaussian  noise 
(AWGN). 

A  secondary  research  objective  wais  to  design  computer  simulation  models  for  each  of  the  code 
tracking  loops  being  evaluated;  such  models  will  be  useful  for  any  follow-on  research. 

1.5  Assumptions 

For  this  thesis,  the  following  assumptions  were  made: 

1.  The  received  GPS  signal  consists  of  a  direct-path  component  and  one  reflected  signal  (i.e., 
two-path  case). 

2.  The  code  tracking  loop  is  already  tracking  the  received  code  prior  to  the  introduction  of  the 
multipath  signal  (i.e.,  the  acquisition  phase  is  complete). 

3.  The  received  GPS  signal  has  been  down-converted  from  the  GPS  carrier  frequency  to  the 
intermediate  frequency  (IF). 

4.  The  received  GPS  signal  corresponds  to  only  one  GPS  satellite  (i.e.,  one  satellite  per  receiver 
channel). 
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5.  Multipath  mitigation  is  achieved  through  signal  processing  alone;  no  special  antennas  or 
spatial  processing  techniques  are  used. 

1,6  Approach 

This  thesis  presents  results  based  on  the  theoretical  analysis  and  computer  simulation  of  the 
NCDLL  and  MRDLL  code  tracking  loop  designs.  Analyses  and  simulations  were  performed  for  a 
variety  of  GPS  direct-path,  multipath,  and  additive  white  Gaussian  noise  (AWGN)  environments. 

All  computer  simulations  are  described  in  Chapter  IV.  The  majority  of  the  code  loop  sim-  * 
ulation  was  performed  using  the  Signal  Processing  Workstation  (SPW®)  software  from  Comdisco 
Systems,  Inc.  of  Foster  City,  California.  Some  preliminary  simulation  required  using  the  Matlab® 
computational  software  from  The  Math  Works,  Inc.  of  Natick,  Massachusetts.  Simulations  were  per¬ 
formed  on  the  Sun  Workstations  provided  at  the  Air  Force  Institute  of  Technology  (AFIT). 
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II.  NCDLL  Analysis 


2. 1  Overview 

This  chapter  presents  an  analysis  of  the  non-coherent  delay-lock  loop  (NCDLL).  The  first 
section  investigates  the  operation  of  the  NCDLL  when  no  multipath  is  present;  this  investigation 
leads  to  the  development  of  the  NCDLL  linear  model.  The  second  section  examines  the  effects  of 
multipath  on  the  NCDLL  with  an  emphasis  on  code  phase  tracking  error  performance. 

A  block  diagram  of  the  NCDLL  is  shown  in  Figure  6.  In  the  analysis  that  follows,  the  following 
assumptions  are  made. 


1.  The  spreading  code  self-noise  is  negligible.  Code  self-noise  is  defined  as  the  time  varying 
component  of  the  NCDLL  discriminator  output,  e  (t,  6).  Generally,  almost  all  of  the  self-noise 
power  is  at  frequencies  which  are  well  outside  the  NCDLL  loop  bandwidth  (5). 

2.  System  processing  gain  ,  G,  is  high.  Processing  gain  is  defined  as 


(8) 


which  is  the  ratio  between  the  GPS  spreading  code  chip  rate,  Wc=  1/%,  and  the  data 
modulation  bit  rate,  Wd-  For  all  spread-spectrum  systems  like  GPS,  Wc  >  Wd\  the  processing 
gain  for  the  C/A  code,  for  example,  is  calculated  as 


Gc/a  — 


(1/Tc) 


c/a  _ 


Wd 


1.023  X  10®  Hz 
50  Hz 


2  X  10\ 


(9) 


and  it  can  be  seen  that  the  processing  gain  for  the  P-code  is  even  greater. 

3.  The  bandpass  filters  (BPFs)  in  Figure  6  are  assumed  to  be  ideal  ‘brick-wall’  filters  centered 
at  IF  with  a  single-sided  noise  bandwidth  of  Bif  Wc  Hz.  This  assumption  gives  a  filter 
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on-time 
code  replica 

C(t-to) 


late 


Figure  6.  Non-coherent  delay  lock  loop  block  diagram, 
transfer  function,  Hipif),  with  magnitude 


\HiFif)\  =  S 


1  for  1/  ±  /cl  <  I  2 
0  elsewhere 


(10) 


Similarly,  the  lowpass  filters  (LPFs)  are  assumed  to  be  ideal  with  two-sided  noise  bandwidth, 
Blpf  ^  Wc  Hz  and  transfer  function 


\HLPF{f)\  = 


1  for  I/I  <  1^1 
0  elsewhere 


4,  The  NCDLL  early-late  chip  offset  has  a  value  A=  1.0.  This  is  the  typical  value  used  in  the 
NCDLL. 

5.  The  BPSK  data  modulation  on  the  received  signal  is  ignored.  This  assumption  is  made 
because  the  NCDLL’s  squaring  circuits  remove  the  effects  of  phase  modulation.  Also,  it  can 
be  shown  that  removing  the  data  modulation  results  in  the  maximum  possible  noise  power  at 
the  NCDLL  phase  discriminator  output  (5),  Therefore,  ignoring  data  modulation  provides  a 
worst-case  scenario  for  the  loop  noise  analysis. 
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6.  The  period,  N,  of  the  spreading  code  is  large.  Therefore,  the  code  autocorrelation  function, 
Rc{il),  is  well-approximated  by  Equation  5  in  Chapter  I. 

7.  Unless  otherwise  stated,  code  Doppler  shift  is  ignored. 

2.2  NCDLL  analysis  with  no  multipath 

2.2.1  The  received  signal.  The  received  direct-path  GPS  signal,  rdp{t),  is  a  modulated 
carrier  in  bandlimited  AWGN;  the  received  signal  (neglecting  data  modulation)  is  given  as 

rdp{t)  =  VWaocit  -  To)  cos  {2Tr fct  +  So)  +  n{t)  (11) 

where  P  is  the  transmitted  signal  power  in  Watts,  uo  is  the  direct-path  attenuation  coefficient,  c(t) 
is  the  DS/SS  BPSK  spreading  code,  tq  is  the  direct-path  signal  propagation  delay  in  seconds,  Wc 
is  the  received  carrier  frequency  in  rad/sec,  and  do  is  the  received  carrier  phase  in  radians  given  by 

^0  =  -WcTo. 

The  noise,  n(t),  is  assumed  to  be  band-limited  zero-mean  AWGN  with  a  two-sided  power 
spectral  density  (PSD)  of  No/2  W/Hz,  and  is  represented  as 

n{t)  =  V2n/(t)  cos  Wet  -  \/2nQ{t)  sin  w^t  .  (12) 

Using  the  representation  in  Equation  12,  the  in-phase  and  quadrature  components  of  the  noise, 
(nj(t)  and  nQ{t),  respectively)  are  independent  zero-mean  lowpass  white  Gaussian  noise  processes 
each  having  a  two-sided  PSD  of  No/2  W/Hz. 

2.2.2  The  early-late  correlator  outputs.  Next,  consider  the  outputs  of  the  early/late 
correlators,  yi{t)  and  j/2(t).  Under  our  previously  stated  assumptions,  it  can  be  shown  that  the 
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signal  components  of  interest  are  (5) 


yi{t)  =  V2PaoRc  ^  ^ Oq) nyi{t)  (13) 

y2{t)  =  V^aoRc  —  2^  +  ^o)  +  '^y2{t) 


where  <5  is  the  normalized  code  phase  estimation  error  defined  cis 


ro(i)  -ro(t) 


and  ro{t)  is  the  NCDLL’s  estimate  of  the  direct-path  propagation  delay.  From  this  point  forward, 
the  implicit  time  dependence  of  S  will  be  dropped  to  simplify  notation. 

The  correlator  output  noises,  nyi{t)  and  ny2{t)^  are  given  by 


nyi{t)  =  c  (i  -  To  +  ^Tc)  n(i)  (15) 

ny2{t)  =  c  (i  -  To  -  ^Tc)  n(t) 


where  n{t)  is  defined  in  Equation  12.  Both  nyi{t)  and  ny2{t)  are  Gaussian  since  n{t)  is  zero-mean 
Gaussian  and  c{t)  equal  -|-1  or  -1  with  equal  probability;  however,  these  correlator  output  noises 
are  obviously  not  independent  (5). 

2.2.3  Discriminator  output.  To  determine  the  NCDLL  discriminator  output,  e  {t,6),  first 
consider  the  signal  component  (neglecting  noise).  This  signal  component  is  obtained  by  squaring 
the  signal  components  of  the  correlator  outputs,  taking  their  difference,  and  taking  the  lowpass 
component  to  give 

eit,6),,y=[yl{t)-ymLP  =  (16) 
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NCDLL  S-curve 


Figure  7.  NCDLL  S-curve. 

where  S  {S)  is  referred  to  as  the  NCDLL  ‘S-curve’  and  is  defined  as 

+  .  (17) 

This  curve  is  plotted  in  Figure  7,  Note  that  there  is  a  linear  region  around  6  =  0.  This  region  is 
always  chosen  as  the  operating  region  of  the  NCDLL,  and  the  loop  will  tend  to  operate  at  the  point 
where  S{S)  =  0  and  the  slope  is  positive  (5).  Figure  7  shows  that  this  operating  point  corresponds 
to  zero  steady-state  code  phase  tracking  error  (i.e.,  633  =  0). 

Next,  consider  the  noise  component,  (i)?  fhe  discriminator  output.  Because  the  loop 
filter  following  the  discriminator  has  a  bandwidth  much  smaller  than  that  of  (t),  only  the  lowpass 
component  is  of  interest.  This  lowpass  noise  component  is  assumed  to  be  AWGN  and  can  be  shown 
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to  have  a  two-sided  PSD,  GncC/)^  given  by  (5) 


Gneif)  : 


Gne{f)\f=o  =  2NqBif  +  2PalNQ  Rl  “  2)  2) 


which  is  defined  over  the  loop  bandwidth  and  is  in  units  of  W/Hz.  The  combined  signal-plus-noise 
expression  for  the  NCDLL  discriminator  output  is 


€  {t,  6)  -  Pals  {S)  +  Tie  (t)  . 


(19) 


2.2.4  NCDLL  linear  model.  The  analysis  in  the  preceding  paragraphs  allows  us  to  form 
a  linear  model  for  the  NCDLL.  This  linear  model  is  valid  for  small  tracking  errors  {S  «  0)  and  is 
formed  by  first  considering  the  operation  of  the  volt  age- controlled  clock  (VCC)  described  by 

I^=gJ\(X,6)*f{X)dX  (20) 

Jo 

where  *  denotes  convolution,  Qc  is  the  VCC  constant  in  Hz/volt,  and  f  (t)  is  the  impulse  response 
of  the  NCDLL  loop  filter.  For  small  tracking  errors  and  large  code  period,  JV,  the  linear  region  of 
the  S-curve  can  be  expressed  via  Equations  17  and  5  as 


(21) 


Therefore,  applying  Equation  21  to  Equation  19,  Equation  20  can  be  rewritten  as 


MO 

Tc 


^ncQc 


C  r  / 

Jo  L  \ 


To  (A)  -  Tb  (A)  A  rie  (A) 


Kn 


fiX)dX 


(22) 


where  Knc  —  2PaQ. 
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Equation  22  represents  a  NCDLL  linear  equivalent  circuit  and  is  shown  in  Figure  8  with 


corresponding  Laplace  transform  shown  in  Figure  9. 


2,2.5  Loop  dynamics  and  tracking  performance.  Examination  of  Figure  9  allows  us  to 
write  the  NCDLL  closed-loop  transfer  function  as 


jj(  \  A  ^q(^)  _  Knc9cF{s) 

n{s)  ~  S  +  KncgcFis) 


(23) 
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In  this  equation,  F{$)  is  the  s-domain  transfer  function  of  the  loop  filter,  which  is  typically  chosen 
to  be  an  active  lead-lag  filler  of  the  form 


F{s) 


1  -hr2S 

TiS 


(24) 


The  presence  of  the  single  power  of  s  in  the  denominator  of  24  allows  the  loop  filter  to  act  as  an 
integrator.  In  this  case,  the  NCDLL  is  called  a  type  2  loop  because  it  contains  two  integrators  (the 
VCC  being  the  second  integrator).  The  advantage  of  a  type  2  loop  is  that  it  can  track  a  signal  in  the 
presence  of  Doppler  shift  (i.e.,  the  introduction  of  a  phase  ramp)  without  incurring  a  steady-state 
phase  error  (2). 

The  closed-loop  transfer  function  can  now  be  rewritten  as  (2) 


H{s) 


2C,WnS  +  UJn 

-h  4- 


(25) 


where  Wn  is  the  loop  natural  frequency  in  rad/sec  given  as 


/ 


KficQc 

n 


(26) 


and  C  is  the  loop  damping  factor  given  as 


/•  '^2 

<:  = 


(27) 


The  single-sided  noise  equivalent  bandwidth  of  the  loop,  Bl  (in  Hz),  can  be  shown  to  be 


Bl  = 


1  t°° 

-j  JH{j2wf)fdf 

i‘^n(c  +  ^)  . 


(28) 
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The  natural  frequency  and  damping  ratio  determine  tracking  performance  of  the  NCDLL  in  the 
presence  of  signal  dynamics  such  as  code  Doppler  shifts;  the  higher  the  natural  frequency,  the 
faster  the  dynamic  response.  The  noise  bandwidth  determines  the  amount  of  noise  power  at  the 
loop  input;  the  smaller  the  bandwidth,  the  greater  the  loop  signal-to-noise  ratio  (SNR).  Equations 
26  to  28  demonstrate  the  fundamental  tradeoff  in  tracking  loop  design:  better  dynamic  response  is 
achieved  only  at  the  expense  of  increased  noise  sensitivity. 

Another  important  NCDLL  performance  parameter  is  the  variance  of  the  steady-state  track¬ 
ing  error,  aj,  also  called  tracking  jitter.  The  power  spectrum  of  the  tracking  jitter,  Gs{f),  is  given 

by  (2) 

Gsif)  =  \Hij2iTf)fG„Af)  (29) 

where  <5„j'(/)  =  Gne{f)l K'^c  ^^e  noise  PSD  at  the  input  of  the  NCDLL  linear  model  in  Figure 
8.  This  noise  PSD  is  approximately  flat  over  the  relatively  small  loop  bandwidth;  therefore,  the 
tracking  jitter  can  be  written  as 


/OO 

Gs{f)df 

-OO 

/OO 

G^Af)m3‘i^f)?df 

-OO 


gn.(0) 

Klc 


Bl 


(30) 


where  Gn€(0)  is  given  in  Equation  18. 


2.S  NCDLL  operation  in  the  presence  of  multipath 

2.3.1  Additional  assumptions.  This  section  investigates  the  operation  of  the  NCDLL  in 
the  presence  of  multipath  interference.  Particular  emphasis  is  placed  on  the  effects  of  the  multipath 
signal  on  NCDLL  steady-state  tracking  performance.  In  addition  to  the  assumptions  stated  at  the 
beginning  of  the  chapter,  two  additional  assumptions  are  necessary: 
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1.  The  purpose  of  this  analysis  is  to  determine  the  effects  of  multipath  interference  on  NCDLL 
operation.  Therefore,  in  order  to  isolate  multipath  effects,  received  signal  noise  is  neglected. 


2.  The  received  multipath  signal  consists  of  a  direct-path  component  with  a  single  reflection. 

2.3.2  The  received  multipath  signal.  The  received  multipath  signal  (neglecting  data  mod¬ 
ulation)  is  given  by 

1 

Tmpit)  =  ajC  {t  -  To  -  aiTc)  COS  {(Vet Oi) n{t)  (31) 

*=o 

where  ai  is  an  attenuation  coefficient  and  Oi  is  the  received  signal  carrier  phase  in  radians  such 
that  Oi  =  -oJe  (to  +  aiTc).  The  parameter,  ai  ,  is  the  delay  coefficient  which  is  defined  such  that 
ao  =  0  and  ai  is  any  real  number  greater  than  zero;  from  this  point  forward,  let  ai  =  a.  We  are 
concerned  only  with  0  <  a  <  1.5  because  relative  multipath  delays  outside  this  region  have  little 
effect  on  GPS  code  phase  tracking  (11).  As  before,  the  received  noise,  n{t),  is  modeled  as  bandpass, 
zero-mean,  AWGN  according  to  Equation  12. 

2.3.3  NCDLL  discriminator  output  in  the  presence  of  multipath.  Neglecting  noise  and 
applying  the  signal  component  of  Equation  31  to  the  NCDLL  input,  the  correlator  outputs  now 
become 


yi{t)  =  ^aoRc  ^ cos  {ufei ^ a)  ^  ^  ^  +  Oi)^ 

y2{t)  =  \/W  |ao/2c  cos  {coct  -h  ^o)  +  «i-^c  ^  “  2)  * 

Therefore,  the  discriminator  output  is 

=  [2/2(0 -2/? (0]lp 

=  P  {ao5  (S)  +  alS{6  +  a)  +  2aoai  cos{au)cTc)  Scrossi^.Oi)} 
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where  Scrossi^,^)  results  from  the  presence  of  cross-terms  and  is  given  by 


Scrosai^,  «)  =  iZc  ^<5  -  -  j  iic  ”  2  )  “  2  )  ^  2  )  ' 

2,S,4  Effect  of  multipath  on  NCDLL  tracking  performance.  Equation  32  is  illustrated 
in  Figure  10.  Note  that  the  NCDLL  discriminator  output  no  longer  consists  of  a  single  S-curve 
centered  about  6  =  0.  Instead,  the  output  is  now  the  summation  of  three  different  S-curves, 
each  attenuated  by  a  different  amount  and  shifted  by  a  different  amount  along  the  6-axis.  The 
net  result,  shown  in  Figure  11,  is  a  tracking  point  that  is  no  longer  at  6  =  0;  this  results  in  a 
steady-state  code  phase  tracking  error  that  depends  on  the  signal  parameters,  ao,  and 

a.  Figure  12  shows  the  predicted  steady-state  tracking  error  for  different  values  of  a  when  ao  =  1, 
ai  =  0.5,  and  the  multipath  component  is  in  carrier  phase  with  the  direct-path  component  (i.e., 
0Q  ^  0i  =  a2irfcTc  =  2n7r  for  n  an  integer).  This  plot  was  obtained  by  setting  Equation  32  equal 
to  zero  for  a  given  a  and  solving  for  6. 

Let’s  examine  the  effect  of  the  phase  difference,  Oq  —0i  =  a2irfcTc^  on  the  NCDLL  multipath 
tracking  error.  It  can  be  seen  that  cos  {a27rfcTc)  in  the  cross-term  of  Equation  32  is  periodic  in 
a  with  period  n/fcTc  (n  an  integer).  In  (12),  Van  Nee  states  that,  for  a  given  set  of  attenuation 
coefficients,  oq  and  ai,  NCDLL  code  phase  errors  are  maximum  when  the  reflected  signal  is  in 
carrier  phase  or  180  degrees  out  of  phase  with  the  direct-path  signal.  To  demonstrate  this,  Figure 
13  plots  the  steady-state  tracking  error  and  the  phase  difference  versus  a  over  approximately  one 
period  of  cos(6>o  -  ^i)  for  ao  =  1,  ai  =  0.5,  and  fcTc  =  10;  note  that  the  peaks  corresponding  to 
maximum  phase  error  occur  when  Oq  —  0i  =  0  and 
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0.25 1 


Predicted  NCDLL  multipath  tracking  error 


Figure  12. 


Figure  13. 


Predicted  NCDLL  multipath  tracking  error  for  ag  =  1,  ai  =  0,5,  and  Oq  -  0i  =  2Ti7r 
(for  n  an  integer). 


Multipath  tracking  error  and  carrier  phase  difference 


Relationship  between  NCDLL  steady-state  tracking  error  and  multipath  carrier  phase 
difference  for  ao  =  1,  ai  =0.5,  and  fcTc  =  10. 
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III.  MRDLL  analysis 


3.1  Overview 

This  chapter  investigates  the  operation  of  the  modified  RAKE  delay-lock  loop  (MRDLL). 
The  first  section  describes  the  received  GPS  signal  and  its  conversion  to  baseband  as  it  enters  the 
MRDLL.  The  last  three  sections  describe  the  operation  of  the  MRDLL’s  three  main  components: 
the  multiple-correlator  tracking  loop  (MCTL),  the  adaptive  loop  controller  (ALC),  and  the  multiple- 
correlator  estimation  unit  (MCEU).  In  the  analysis  that  follows,  the  assumptions  presented  at  the 
beginning  of  Chapter  II  still  apply. 

3.2  Received  signal  model 

The  MRDLL  block  diagram  is  shown  in  Figure  14.  As  in  Chapter  11,  the  received  GPS 
multipath  signal  is  modeled  as  having  a  direct-path  component  and  a  single  reflection  according  to 
Equation  31.  After  conversion  to  baseband  as  depicted  in  Figure  14,  the  signal  into  the  MRDLL 
becomes 

r(t)  =  xoc(t  -  To)  -I-  Xic{t  -To  -  aTc)  +n'{t)  (34) 

where  the  coefficients,  xq  and  xi,  are  defined  as  {O3  is  the  locally-generated  carrier  phase) 

xo  =  \/Waocos{Oo  -  63)  (35) 

xi  =  V^ai  cos  {$1  -  O3)  , 

and  n'{t)  is  a  lowpass,  zero-mean,  AWGN  process  given  as  (see  Equation  12) 

n'(t)  =  v/2[n/(t)cos(^3)  -I- re(3(f)  sin  (^3)]  (36) 

which  has  a  two-sided  PSD  of  Nq  W/Hz  over  the  LPF  bandwidth,  B^pp  Hz. 


26 


MCTL 

Figure  14.  Modified  RAKE  delay  lock  loop  (MRDLL)  block  diagram. 

As  shown  in  Figure  15,  the  local  carrier  is  generated  using  a  code-correlator  and  phase-locked 
loop  (PLL)  in  combination  (4).  To  derive  an  expression  for  the  local  carrier  phase,  ^3,  we  assume 
negligible  code  phase  tracking  error  (i.e.,  fo  =  tq)  and  neglect  the  effects  of  noise  on  the  PLL  (i.e., 
O2  =  ^0  when  no  multipath  is  present).  The  signal,  into  the  PLL  can  be  written  as 

r\t)  =  VWao  cos{wct  +  ^0)  +  V^aiRc  (a)  cos  {coct  +  ^1)  (37) 

where  Oq  -  0i  =  2'KfcOtTc^  Using  trigonometry,  Equation  37  can  be  rewritten  as 

r\t)  =  C  cos  {ujct  +  ^0  +  0e)  (38) 

where  the  phase  error,  due  to  the  reflected  signal  is  a  function  of  o.  and  flo/^i  ^nd  is  described 

by 

_i  [  -i?e(Q!)sin(<9o -6>i)  1  ^  r  -Re(Q=)sin(27r/eaTo)  1  . 

^  ao/ai Rc{oi)  cos  {Oq  —  6 1)  _aQ / ai Rc  (a)  cos  {2'ir f cOiTc) 
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2  cos  (cOcrH-03 ) 
(to  baseband  mixer) 

c(/-To) 

Figure  15.  Typical  carrier  phase  recovery  scheme  for  GPS. 


Figure  16.  Local  carrier  phase  error,  <^e,  for  a^/ai  =  2  and  fcTc  =  10. 

Therefore,  the  phase  tracked  by  the  PLL  (and  hence,  the  phase  of  the  locally-generated  carrier)  is 

^3  ~  ^0  +  0e- 

Figure  16  plots  (j>e  versus  a  for  fcTc  —  10  and  a^/a^  =  2;  note  that  the  local  carrier  is  in  phcise 
with  the  direct-path  carrier  (i.e.,  (f>e  =  0)  whenever  a  takes  on  values  for  which  the  reflected  signal 
is  in  phase  or  180  degrees  out  of  phase  with  the  direct  path  component.  Also,  note  from  Figure 
16  that  the  phase  error  envelope  follows  the  code  autocorrelation  function,  and  the  local 

carrier  phase  error  is  zero  for  all  a  >  1;  this  is  due  to  the  fact  that  the  code  correlation  preceding 
the  PLL  effectively  removes  the  reflected  signal  component  (i.e.,  i2c(<^)  =  0  for  Oi  >  1), 
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To  see  the  effect  oi4>e  on  the  signal  parameters,  lo  and  xi,  we  can  use  trigonometric  identities 


% 


and  Equation  35  to  write  the  ratio,  xq/xi  as 


xq  ^  ao 

Xi  dl 


(40) 


where 

_ I - 1  .  (41) 

^  [cos  (ucOtTc)  -  sin  (wcaTc)  tan  0*  J 

Note  that  for  zero  phase  error,  4>e-,  the  ratio  in  Equation  40  reduces  to  ±ao/ai  depending  on  whether 
the  reflected  signal  is  in-phase  (i.e,  2'irfcaTc  =  2n7r)  or  out  of  phase  (i.e.,  2'KfcO.Tc  =  (2n  -1-  l)ir) 
with  the  direct-path  signal  for  n  an  integer. 


S.3  MCTL  operation 

3.3.1  Overview.  The  MCTL  is  responsible  for  tracking  the  direct-path  code  phase  of 
the  incoming  multipath  DS/SS  code.  Like  the  NCDLL,  the  MCTL  uses  early-late  code  correlation 
(without  the  squaring  operation)  to  perform  tracking.  However,  unlike  the  NCDLL,  the  MCTL 
relies  on  signal  parameter  estimates  to  accurately  track  the  code  phase  in  the  presence  of  multipath. 
Therefore,  the  MCEU  has  been  designed  to  provide  the  MCTL  with  maximum-likelihood  (ML) 
estimates  of  the  multipath  delay  coefficient,  a,  and  of  the  signal  parameters  xo  and  xi-  As  will 
be  seen  later,  these  estimates  enable  the  MCTL  to  remove  the  tracking  error  introduced  by  the 
reflected  signal. 

As  shown  in  Figure  14,  the  MCTL  consists  of  Af  correlator  arms.  In  the  kth  arm  of  the 
MCTL,  the  received  signal  is  correlated  with  shifted  early  and  late  replicas  of  the  DS/SS  code.  For 
an  early-late  spacing  of  one  chip  length,  these  shifted  replicas  are  defined  as 

CA  (t  -  TO  -  /3fcTo)  ^  c  (f  -  fo  -  HkTo  -  Tc/2)  -  c  (t  -  fo  -  l3kTc  +  Ta/2)  (42) 
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where  fik  is  the  delay  spacing  coefficient  of  the  kth  loop  arm  and  tq  is  the  MCTL’s  estimate  of  the 
direct-path  delay.  The  arm  corresponding  to  =  0  is  referred  to  as  the  direct-path  arm  (/3o  =  0) 
and  the  remaining  arms  as  multipath  arms.  During  tracking,  the  MCTL  enables  both  the  direct- 
path  arm  and  the  multipath  arm  corresponding  to  =  3  (where  a  is  the  MCEU’s  estimate  of  the 
multipath  delay  coefficient). 

In  the  analysis  that  follows,  assume  that  the  the  MCEU  is  providing  perfect  estimates  to  the 
MCTL  and  that  the  multipath  arm  corresponding  to  =  1  is  enabled.  In  other  words, 


/3i  =  a  =  O!,  xq  =  xq,  and  xi  =  xi  .  (43) 

3.3.2  Early-late  correlator  outputs.  The  baseband  signal  undergoes  early-late  correla¬ 
tion  in  each  of  the  enabled  arms.  Under  our  assumption  of  exact  MCEU  estimates,  the  lowpass 
components  of  the  correlator  outputs  are 

zo(t,^)  =  xqD  {6)  +  X\D  {6 a) (44) 

zi(t,  6)  =  xqD  {6  -  a) xiD  {6) 

where  the  MCTL  D-curve.,  -^(^)i  i®  defined  as 

Z)(5)4/2e(<5-i)  -i2e(^+^)  .  (45) 

Note  that  this  curve  is  almost  identical  (except  for  the  squaring  operation)  to  the  NCDLL  S-curve 
of  Equation  17.  The  noise  components,  770  (0  and  77i(t),  are  described  by 
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Tfl{t)  =  n/(i)cA  it-%  -  PkTc)  *  hipf{t) 

Vki^)  —  (t-To  -  PkTc)  *  hipf{t) 

where  *  denotes  convolution  and  hipf(t)  is  the  impulse  response  of  the  LPFs.  It  can  be  shown  that, 
for  code  period  iV  »  1,  the  two-sided  PSD,  G,,' ,  of  r}l{t)  and  is  approximately  (5) 


Gr,',{f)  =  No 


(47) 


in  units  of  W/Hz.  Therefore,  the  output  noise  of  the  fcth  LPF  given  by  Equation  46  has  a  two-sided 


PSD  of 


2iVo  I/I  < 


Gr,Af)=< 


(48) 


0  elsewhere 


3.3.3  Discriminator  output.  After  early-late  correlation,  the  signal  enters  the  gain-phase 
correlators.  Here,  the  signal  is  mixed  with  the  appropriate  MCEU  signal  parameter  estimate,  xq  or 
xi,  as  shown  in  Figure  14.  Assuming  perfect  estimation,  the  outputs  of  the  gain-phase  correlators 
are 


yo{t,6)  =  xlD  (6)  +  xqXiD  {6  +  a)  +  xoVoit)  (49) 

yi(t,6)  =  xqXiD  (6  -  a)  +  xlD  {S)  +  xir]i{t)  . 


These  correlator  outputs  are  summed  to  produce  the  MRDLL  discriminator  output,  e{t),  given  by 


e(t,(5)  =  yo{t,6)  +  yi{t,S)  =  5 (6)  +  ne{t)  ; 


(50) 


the  MRDLL  S-curve  is  given  as 

S  (6)  =  (a:o  +  xf)  D  (5)  +  xqXi  \D  {6  a)  +  D  {6  -  a)]  ,  (51) 

and  the  discriminator  output  noise,  ne(t),  is  given  by 

ne(t)  =  XoTJoit)  +  XtT)i{t)  .  (52) 


Using  the  result  of  Equation  48,  the  discriminator  output  noise  PSD  becomes 


Gneif)  = 


{xl+x\)  (2iVo) 

0 


l/l<  ^ 

elsewhere 


(53) 


3.3.4  MRDLL  steady-state  tracking  error  due  to  multipath.  To  determine  MRDLL  steady- 
state  tracking  error  due  to  multipath  (neglecting  noise  effects),  consider  the  representative  MCTL 
S-curve  shown  in  Figure  17.  As  was  the  case  for  the  NCDLL  in  the  presence  of  multipath,  the  MCTL 
S-curve  is  the  summation  of  three  separate  curves  multiplied  and  shifted  by  different  amounts.  One 
D-curve  is  centered  at  6  =  0;  the  other  two  D-curves  are  shifted  by  a  in  opposite  directions  along  the 
^-axis.  However,  because  of  the  way  we  chose  to  implement  the  gain-phase  correlation  in  Equation 
49,  the  two  shifted  D-curves  are  each  multiplied  by  the  same  amount.  As  shown  in  Figure  18, 
this  implementation  ensures  that  the  loop  tracking  point  remains  at  6  =  0.  Therefore,  unlike 
the  NCDLL,  the  MRDLL  exhibits  zero  steady-state  code  phase  tracking  error  in  the  presence  of 
multipath  interference  (neglecting  noise  effects). 
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3.3.5  MCTL  linear  model.  The  MCTL  linear  model  for  small  tracking  errors  ((5  «  0)  is 
formed  by  first  considering  the  operation  of  the  VCC,  described  by 


Tc 


=  gc  f  e{\,S)*fiX)dX  , 
Jo 


(54) 


where  *  denotes  convolution,  gc  is  the  VCC  constant  in  Hz/ volt,  and  /  (t)  is  the  impulse  response 
of  the  MCTL  loop  filter  (contained  in  the  ALC).  For  small  tracking  errors,  the  linear  region  of  the 
S-curve  can  be  expressed  as 

S{6)  =  AS  =  a(^'^^^  ,  (55) 

where  A  is  defined  as  the  slope  of  the  linear  operating  region  for  small  6.  As  will  be  seen  in  Section 
3.4,  this  slope  depends  on  the  signal  parameters,  xq,  Xi,  and  a.  Applying  Equation  55  to  Equation 
50,  the  VCC  equation  can  be  rewritten  as 


To(0 


+ 


Ue  (A) 


*  /  (A)  dX 


(56) 


Note  that  Equation  56  is  almost  identical  to  Equation  22  in  Chapter  II;  therefore,  the  linear 
equivalent  circuit  and  the  Laplace  transform  linear  model  for  the  MCTL  can  be  represented  by 
Figures  8  and  9,  respectively,  with  Knc  replaced  by  A,  and  n^lt)  replaced  by  ne(t). 


3.3.6  MCTL  dynamics  and  tracking  performance.  The  MCTL  closed-loop  transfer  func¬ 


tion  can  be  written  as 

ro(s)  s  +  AgcF{s) 


(57) 


where  F{s)  is  the  transfer  function  of  the  loop  filter.  The  MCTL  is  chosen  to  be  a  type  2  loop 
with  the  active  lead-lag  loop  filter  transfer  function  described  by  Equation  24;  the  loop  natural 
frequency  in  rad/sec  is  given  as 


Wn  = 


(58) 
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and  the  damping  ratio  and  noise  bandwidth  are  described  by  Equations  27  and  28. 

The  power  spectrum  of  the  tracking  jitter  is  given  by  (2) 

Gsif)  =  \Hij27tf)fG„Af)  (59) 

where  G„j'(/)  =  Gne{f)/A'^  is  the  noise  PSD  at  the  input  of  the  MCTL  linear  model.  This  noise 
PSD  is  approximately  flat  over  the  relatively  small  loop  bandwidth;  therefore,  the  tracking  jitter 
can  be  approximated  as 

=  r  Gs{f)df  =  r  G„,>  if)  \Hij27rf)f  df  =  (60) 

J —oo  J —  oo 

i{xl  +  xl)No 

=  - ^2 - 

where  Gne(O)  was  determined  from  Equation  53. 

3.4  Adaptive  loop  controller  (ALC)  operation 

The  block  diagram  for  the  ALC  is  shown  in  Figure  19.  The  ALC  (assuming  a  second  order 
loop)  consists  of  a  loop  filter  defined  by 


fm  =  (61) 

preceded  by  a  variable  gain,  Ka=  I/ti;  this  effectively  gives  a  loop  filter  trajisfer  function  identical 
to  Equation  24  (i.e.,  F{s)  =  As  will  be  seen,  the  variable  gain  allows  the  loop  designer 

to  fix  the  dynamic  performance  of  the  MCTL  about  the  linear  operating  region  (i.e.,  assuming 
small  6).  The  ALC  feature  is  unique  to  the  MRDLL  and  was  not  present  in  the  RDLL  design  in 

(7). 
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Figure  19.  Adaptive  loop  controller  (ALC)  block  diagram. 

Adaptation  to  the  multipath  environment  enables  consistent  dynamic  code  tracking  perfor¬ 
mance  in  the  presence  of  code  Doppler  shift.  Such  a  Doppler  shift  is  naturally  created  by  satellite 
motion  and  can  be  intensified  by  aircraft  maneuvers.  The  need  for  adaptive  control  is  embodied  in 
Equation  58.  This  Equation  shows  that  the  MCTL  natural  frequency  (and,  therefore,  the  damping 
ratio  and  noise  bandwidth)  is  dependent  upon  the  value  of  the  S-curve  slope,  A,  which  in  turn 
depends  on  the  signal  parameters,  xq,  xi,  and  a;  the  MCTL  linear  operating  region  is  defined  in 
Table  2. 


Table  2.  MCTL  linear  operating  region. 


Gain  =  (A/2) 

Multipath  Delay  Range 

Tracking  Error  Range 

XQ  +  x'f  +  2xoXi 

0  <  a  <  0.5 

—0.5  +  a  <  <5  <  0.5  —  a 

Xq  +  -h  0,5xoXi 

a  =  0.5 

-0.5  <  5  <  0.5 

xl  +  x'(-  XoXi 

0.5  <  a  <  1 

1  0.5  —  a  <  ^  <  -0.5  +  a 

xl-\-x{-  IflXl 

1  <  Q!  <  1.5 

-1.5  + a  <  6  <  1.5-q 

Xq  -f  x^  -  0.5xoXi 

a  =  1.5 

-0.5  <S<0.5 

Table  3.  MCTL  affine  tracking  region. 


Slope 

Multipath  Delay  Range 

Tracking  Error  Range 

Xq  +  xf  +  0.5xoXi 

0  <  a  <  0.5 

|5|  >  0.5-a 

Xq  +  Xi  +  0.5xoXi 

0.5  <  a  <  1 

|5|  >  —0.5  +  O' 

Xq  +  x^  -  0.5xoXi 

1  <  a  <  1.5 

|5|  >  1.5-a 
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a  =  0.4 


Figure  20.  MRDLL  S-curve  with  xq  =  1,  =  0.5,  and  a  =  0.4. 

Figures  20  through  24  show  typical  MCTL  S-curves  for  different  values  of  a  with  xq  =  1 
and  xi  =  0.5;  note  the  change  in  slope  in  the  linear  region  about  the  origin.  Also  note  (for 
a  ^  0.5,  1.5)  that  outside  this  linear  region  (into  the  affine  region  defined  in  Table  3),  the  S-curve 
slope  slightly  changes,  but  the  overall  S-curve  for  \S\  <  0.5  still  remains  approximately  linear; 
therefore,  loop  dynamic  performance  will  vary  in  the  affine  tracking  region,  but  only  slightly  from 
that  of  the  linear  region.  Both  the  figures  and  Table  2  demonstrate  a  significant  difference  between 
the  MRDLL  and  the  RDLL.  In  the  RDLL,  the  equation  governing  the  linear  operating  region  does 
not  change  with  integer  chip  delays  and  is  always  defined  over  -0.5  <  (5  <  0.5  (7). 

To  maintain  a  fixed  natural  frequency,  (and  therefore  a  fixed  damping  ratio  and  noise 
bandwidth),  the  ALC  adapts  the  value  of  ri  via  estimates  from  the  MCEU;  this  approach  is  based 
on  the  certainty  equivalence  principle  (1)  wherein  the  controller  is  designed  using  estimates  of 
unknown  parameters.  The  order  of  operation  for  the  ALC  is 

1.  Receive  signal  estimates  (xo,  xi,  and  a)  from  the  MCEU 
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(S)S' 


2.  Use  the  MCEU  estimates  to  compute  an  estimate,  according  to  Table  2 

3.  Use  A  to  adjust  the  gain,  Ka^  according  to 


=  (62) 

n  Age 

which  was  obtained  by  rearranging  Equation  58  and  substituting  A  for  A. 

One  final  design  consideration  needs  to  be  made.  For  a  second-order  loop  to  operate  correctly, 
the  relationship 

Age  >  (63) 

must  be  satisfied  (2).  This  condition  can  be  restated  by  applying  Equation  58  to  give 

Ka  =  —  •  (64) 

n 

Therefore,  the  ALC  must  be  designed  so  that  the  variable  gain,  Ka,  does  not  exceed  some  threshold. 
If  we  represent  this  gain  threshold  as  ATjimax)  then 


Ka  <  KAn^^x  (65) 

must  be  satisfied  for  the  MCTL  to  operate  correctly.  The  value  of  KAm&x  i®  chosen  by  the  loop 
designer  to  optimize  MCTL  performance. 

3.5  MCEU  Analysis 

3.5.1  Overview.  We  have  already  seen  that  if  the  MCEU  provides  perfect  estimates 
to  the  MCTL,  then  we  are  able  to  remove  steady-state  tracking  error  due  to  a  reflected  signal; 
furthermore,  these  estimates  allow  the  ALC  to  maintain  consistent  dynamic  performance  in  the 
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Xi 


T  c(t-  vPm- Jc) 

Figure  25.  Multiple-correlator  estimation  unit  (MCEU)  block  diagram. 

MCTL  (assuming  linear  operation).  Now,  we  will  examine  the  operation  of  the  MCEU  and  the 
theory  behind  it. 

S.5.2  Sampling  the  correlation  function.  The  MCEU  block  diagram  is  shown  in  Figure  25. 
The  MCEU  consists  of  two  components:  a  correlator  bank  and  an  estimator.  The  correlator  bank 
consists  of  M  arms  in  which  the  A:th  arm  is  correlated  with  a  replica  code  given  by  c{t -tq  -  PkTc). 
Assuming  that  loop  tracking  errors  are  negligible  (i.e,,  tq  =  tq),  the  output  of  the  kth  correlator 
into  the  estimator  is  given  by 

Rk{t)  =  xoRciPk)  “I"  xiRc{cx  —  (3k)  4*  Vk{t)  .  (66) 

The  correlator  output  noise,  Vk{t),  assumed  to  be  zero-mean,  lowpass  AWGN,  is  given  by 

Vk{t)  =  n'{t)c{t  -  To  -  f3kTc)  *  hipf{t)  (67) 

where  n'{t)  is  the  received  baseband  signal  noise  defined  in  Equation  36.  This  noise  is  assumed  to 
have  a  two-sided  PSD  of  Nq  W/Hz  over  Bipf  Hz. 
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Examining  the  signal  component  of  the  correlator  output  in  Equation  66;  this  expression 
shows  that,  in  effect,  the  A;th  correlator  arm  ‘samples’  both  the  direct  path  component  and  the 
corresponding  multipath  component  of  the  code  correlation  function,  iZc(^i)»  at  fi  =  The 
estimator  then  uses  the  correlator  bank  outputs  as  discrete-time  measurements  to  estimate  the 
unknown  parameters,  xq,  Xi,  and  a.  This  concept  of  sampling  the  correlation  function  is  presented 
in  (9);  however,  the  analysis  in  (9)  does  not  specify  the  estimation  scheme  employed.  For  our 
estimator,  we  have  chosen  to  use  the  ML  estimation  scheme  presented  in  the  next  section. 


3.5,3  The  MCEU  estimator.  To  begin  the  estimator  analysis,  let  one  sample  of  the 
unknown  signal  parameters  be  represented  in  vector  form,  x,  as 


X  = 


Xi 


(68) 


This  allows  us  to  represent  the  sampled  correlator  outputs,  at  a  given  time  instant,  in  vector  form 
as 

R'  =  [  i?o  Ri  •  *  Rm~i 


—  H(Qi)x  +  V 


(69) 


where  R  is  an  M  x  1  vector  of  the  correlator  outputs,  Rk,  and  H(q:)  is  an  M  x  2  regressor  matrix 
given  by 


U{a)  = 


Rc{M  Rc{oi  —  /3o) 

Rci^i)  Rcia-Pi) 


(70) 


Rc{Pm-i)  Rc{<^  -  Pm-i)  J 

The  M  X  1  noise  vector,  [vo  vi  •  •  •  consists  of  the  correlator  noise  output  at  the  time 

the  measurements,  are  taken  and  is  distributed  as  N  :  [0,  C^],  where  is  the  noise  covariance 
matrix.  To  determine  an  expression  for  Cv,  we  treat  the  spreading  code  as  a  random  binary  process 
and  assume  it  is  independent  of  the  noise;  this  gives  the  continuous-time  cross-correlation  between 
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the  ith  and  jth  correlator  noise  outputs  as 


Rijir)  =  E\vi{t)vj{t  +  T)]  (71) 

=  E  {K(t)c(i  -  fo  -  jiiTc)  *  [n'{t  +  T)c{t  +  r  -  fo  -  ^jTc)  *  hipfit  +  r)]} 

=  E  [c{t  -  fo  -  ATc)c(i  +  T  -  fo  -  ^jTc)\  E  [n"it)n"it  +  r)] 


where 

n"(i)  =  n'(t)  *  hipf{t) 


(72) 


and 

.  (73) 

Because  the  spreading  code  bandwidth  is  nauch  greater  than  bandwidth  of  the  lowpass  noise,  Rn"  (r ) 
is  approximately  constant  over  significant  values  of  Rc  i  the  cross-correlation  can  be 

approximated  by 

Rijir)  «  i?e(r  A/3i,)J?„K0)  =  Rc{r  +  (74) 

where  is  the  variance  of  each  the  correlator  noise  outputs.  This  result  gives  the  discrete-time 
noise  covariance  matrix,  Cv,  as 

(75) 

where  C  is  the  M  X  M  correlation  matrix 

1  Rc{l3o-l3i)  •••  Rc{Po-Pm-i) 

RciPi-lSo)  1  : 

:  Rc{Pm-2  — 

Rc{Pm-1-Po)  •••  Rc{Pm-\  -  Pm~2)  1 
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which  is  symmetric  because  Rc  (A  -  /Jj)  =  Rc  {Pj  -  ^i)  and  Toeplitz  since  Rc  (A+fc  -  ft)  = 

RciPk  -  ft)  for  all  i,  k. 

Assuming  that  a  reflected  signal  is  present  (i.e.,  a  ^  0),  the  maximum-likelihood  (ML) 
estimate  of  x  for  a  given  value  of  a  is  given  by  (6) 

X  =  P(o!)R  (77) 


where 

P(a)  4  [H'^(a)C-iH(a)]  H'^(a)C-^  ;  (78) 

observe  that  knowledge  of  cr^  is  not  needed  to  compute  the  estimate,  x,  of  Equation  77. 

It  is  not  feasible  to  compute  x  per  Equation  77  for  the  continuum  of  a  €  [0,  1.5].  We 
alternatively  choose  to  consider  the  estimates,  Xj,,  generated  by  assuming  that  a  =  ft  for  k  = 

(0,  1,  2,  ...,  M-1).  We  are,  in  effect,  considering  a  flnite  collection  of  uniformly  spaced  possibilities 

for  a.  To  compute  the  estimate,  x*  corresponding  to  the  event  a  =  ft,  the  MCEU  forms  P{^k) 
according  to  Equation  78.  Equation  77  is  then,  for  k  =  (1,  2,  ...,  M  -  1),  used  to  compute  x*.  as 

X*  =  P(ft)R  (fe  5^  0)  .  (79) 


This  computation  is  simplified  by  realizing  that  P(ft)  is  a  2  x  16  matrix  with  non-zero  entries 
only  in  the  first  and  {k  -|-  l)st  columns  (zeros  elsewhere);  furthermore,  the  elements  in  the  (k  -h  l)st 
column  are  simply  the  reverse  of  the  elements  in  the  first  column;  therefore,  the  fcth  estimate  can 


be  rewritten  as 


Xjfc  = 


Pll  P2\ 
P21  Pll 


R, 

Rk 


=  PkRk  (fe^o) 


(80) 


where  Pk  is  a  symmetric  2x2  matrix  with  the  same  first  column  as  P(/3k)  is  a  2  x  1 

column  vector  with  elements  corresponding  to  the  correlator  measurements  Rq  and  Rk> 
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To  compute  the  direct-path  only  estimate  corresponding  to  A;  =  0  (the  event  a  =  0),  we  must 
first  realize  that  H*^(a)C“^H(a)  in  Equation  78  is  singular  when  o:  =  0;  this  means  P{/3o)  =  P(0) 
cannot  be  evaluated  as  in  Equation  79.  Since  we  have  assumed  no  multipath,  we  can  also  assume 
that  =  0  (see  Equation  35);  therefore,  we  only  need  to  estimate  the  direct-path  signal  parameter, 
xq.  This  direct-path  estimate  is  calculated  by  defining  a  constant  16  x  1  vector,  H(^o)?  as 


H(/?o)  ^ 


Rc  (/3o) 

RM 

Rc  (i^M-l) 


and  substituting  this  vector  for  H(a)  in  Equation  78  to  form 


(81) 


P(/3„)  4  [H(;0orC-iH(^o)]  ^H(/3orC-i  ; 


(82) 


this  gives  the  k  =  0th  estimate  as 


Xo  = 

Xq 

_ 

P()0o)R 

Xi 

0 

0 

(83) 


As  for  the  multipath  case,  this  computation  can  be  simplified  by  realizing  that  P(jdo)  is  a  fixed 
1  X  16  vector  with  a  1  in  the  first  element  and  zeros  in  all  the  remaining  elements;  therefore,  the 
direct  path  signal  estimate  is  simply  the  0th  correlator  output  (i.e.,  xq  =  Rq)  because,  as  seen  in 
Equation  66,  Rq  —  xoi?c(0)  and,  from  Equation  5,  Rc{0)  =  1.  In  summary,  the  estimation  vector 
for  the  direct-path  only  case  becomes 


xo 


Rq 

0 


(84) 
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Once  all  estimates  have  been  computed,  the  MCEU  must  determine  which  estimate,  x*,, 
most  likely  gave  rise  to  the  correlator  measurement  matrix,  R;  the  corresponding  /3k  becomes  the 
ML  estimate  for  a  (i.e.,  a  =  /?*,).  To  accomplish  this,  the  MCEU  forms  a  matrix  of  predicted 
observations,  Rk,  for  each  of  the  M  correlator  arms  as 


Rk  =  H(a)xk  =  H(/?k)xk  • 


(85) 


The  prediction  error  vector,  et,  is  then  computed  according  to 

©k —  R.“R-k  * 


(86) 


Consider  the  case  when  the  estimates  equal  the  actual  values;  define  this  event  to  be 


=  {(a  =  a  =  /3fc)  n  (xo  =  xq)  n  (xi  =  xi)}  . 


(87) 


We  want  to  choose  the  most  likely  value  of  ©k  given  that  is  true;  in  other  words,  we  want  to 
maximize  the  likelihood  function,  p  (ek|'®')  • 

Applying  Equations  69  and  85,  and  our  definition  of  ©k  in  Equation  86,  we  get 

©k|^=  R-Rk  =  H(a)x  -  H(3)xk  +  v  =  v  ,  (88) 


because,  under  our  assumed  condition,  H(q!)x  =  H(a)xk.  Therefore,  (©kl^)  is  distributed  as 
N  :  (0,  Cv]  and  this  gives  the  likelihood  function  as  (6) 


p(©k|4')  = 


{27r)^^^  Vdet(Cv) 


exp 


(89) 
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For  a  given  noise  variance,  it  can  be  seen  from  Equation  89  that  the  likelihood  function  is 
maximized  when  eJ^’C^^ejc  is  minimum;  therefore,  the  MCEU  computes  Ek  ~  and 

chooses  the  estimates  corresponding  to  the  minimum  value  of  Ek  (for  A:  =  0,  1,  M  -  1).  These 
estimates  are  sent  to  the  MCTL  as  the  ML  estimates,  x  =  Xjb  and  a  — 
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IV.  Computer  Simulation  Results 


4.1  Overview 

This  chapter  presents  the  results  of  digital  computer  simulations  of  the  NCDLL  and  MRDLL. 
Computer  simulations  were  performed  to  verify  the  predicted  theoretical  results  presented  in  Chap¬ 
ters  I  through  III,  and  to  investigate  actual  loop  performance  in  the  absence  of  simplifying  assump¬ 
tions  such  as  perfect  estimation  and  linear  operation.  Specifically,  four  different  simulations  were 
performed: 

1,  MCEU  estimator  performance  assuming  MCTL  perfect  code  phase  estimation  in  the  presence 
of  AWGN  (preliminary  simulation  using  MATLAB) 

2,  NCDLL  vs  MRDLL  steady-state  tracking  performance  in  the  presence  of  multipath  without 
AWGN 

3.  NCDLL  vs  MRDLL  steady-state  tracking  performance  in  a  direct-path  only  environment 
without  AWGN 

4.  MRDLL  steady-state  tracking  performance  in  the  presence  of  multipath  with  AWGN 

Unless  otherwise  stated,  all  simulations  were  performed  with  the  Signal  Processing  Work¬ 
station  (SPW)  software  package,  and  the  resulting  data  was  analyzed  using  MATLAB,  All  SPW 
computer  models  and  MATLAB  m-files  used  for  simulation  and/or  data  processing  are  contained 
in  Appendices  A  and  B. 

The  remainder  of  this  chapter  is  divided  into  six  sections.  Section  4.2  presents  the  main 
simulation  parameters  and  discusses  how  and  why  the  parameters  were  chosen.  Next,  section  4.3 
describes  the  data  processing  methods  that  were  used.  The  remaining  four  sections  present  the 
objectives,  methods,  and  results  of  each  of  the  four  simulation  scenarios  outlined  above. 
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4-2  Simulation  Parameters 


4.2.1  Simulation  parameters  versus  actual  GPS  parameters.  We  desired  simulation  pa¬ 
rameters  that  reflected  (as  best  as  possible)  actual  GPS  signal  and  tracking  loop  specifications; 
however,  because  SPW  processing  time  was  relatively  slow,  simulation  times  and  data  file  size  be¬ 
came  limiting  factors.  Because  of  this,  we  were  constrained  to  choose  parameters  with  the  intent  of 
demonstrating  general  loop  behavior  and  tendencies  as  opposed  to  attempting  to  simulate  actual 
performance  in  the  GPS  environment;  despite  these  constraints,  simulation  results  still  provided 
valuable  insight  into  the  application  of  the  MRDLL  to  the  GPS  multipath  mitigation  problem.  In 
this  section,  we  look  at  several  important  simulation  parameters,  how  these  parameters  relate  to 
actual  GPS  values,  and  each  parameter’s  potential  effect  on  simulation  results. 

4.2.2  Sampling  frequency.  Three  main  factors  were  considered  in  choosing  the  SPW 
simulation  sampling  frequency,  fs  (in  Hz);  code  phase  timing  error  resolution,  the  number  of 
relative  multipath  delay  samples,  and  simulation  time.  Estimating  the  code  phase  timing  error 
involved  comparing  leading  edges  of  the  loop  on-time  replica  and  the  received  direct-path  code.  The 
relative  time  delay  (or  advance)  between  any  two  leading  edges  had  to  be  estimated  with  adequate 
resolution;  therefore,  a  sampling  frequency  corresponding  to  100  samples/chip  was  chosen.  The 
second  constraint  on  fs  was  that  it  had  to  be  chosen  so  that  the  value  of  the  relative  multipath 
delay,  aTc,  corresponded  to  an  integer  number  of  samples  (i.e.,  aTcfs  had  to  be  a  positive  integer); 
in  our  simulations,  a  took  on  values  equal  to  integer  multiples  of  0.1  between  0  and  1.5.  Taking  all 
these  considerations  into  account,  along  with  simulation  time  and  data  file  size,  led  to  a  choice  of 
fs  =  1000  Hz;  this  corresponded  to  a  chip  rate  of  Wc  =  1/Tc  =  10  Hz  (to  meet  the  100  samples/chip 
requirement). 

4.2.3  Carrier  phase  parameter.  Since  the  NCDLL  code  phase  error  due  to  multipath  is 
maximum  whenever  the  direct-path  signal  is  in  carrier  phase  with  the  reflected  signal  (see  Chaptei 
II),  simulations  were  run  with  =  ^1;  this  allowed  us  to  demonstrate  the  MRDLL’s  ability  to 
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reduce  multipath  tracking  errors  even  in  a  multipath  environment  that  would  otherwise  cause 
significant  tracking  errors  in  the  standard  NCDLL.  Note  that  the  in-phase  condition  also  implies 
that  the  MRDLL  local  carrier  phase  error,  4>ei  described  in  Chapter  III,  is  equal  to  zero;  note, 
however,  this  does  not  imply  a  best-case  scenario  for  the  MRDLL  since  any  local  carrier  phase 
error  is  accounted  for  by  the  signal  estimates,  xq  and  xi  (see  Equation  35). 

4.2.4  Carrier  frequency.  Simulations  were  run  independent  of  choosing  a  carrier  frequency 
parameter,  Uc-  This  was  possible  because  the  complex  envelope  model  described  in  Appendix  A 
is  independent  of  carrier  frequency  except  for  the  carrier  phases,  Oi,  which  depend  on  a  and  Uc- 
Since  we  are  interested  in  relative  multipath  delays,  we  chose  the  direct-path  delay  as  tq  =0,  which 
meant  =  0  and  =  -WcOiTc  =  -2'irfcaTc.  For  our  simulation,  we  assumed  the  product,  fcTc, 
was  large  enough  so  that  afcTc  was  always  an  integer  for  our  chosen  values  of  a  =  O.lfc  and  k  an 
integer  between  0  and  15;  for  comparison,  recall  that  for  GPS,  fcTc  is  a  relatively  large  number 
(e.g.,  fcTc  =  1540  for  C/A  code)  (see  Table  1).  Therefore,  under  our  assumption  of  large  fcTc, 
our  simulated  reflected  signal  carrier  phase,  ^1,  was  always  an  integer  multiple  of  27r.  This  was 
equivalent  to  our  setting  60  =  0i  during  simulation. 

4.2.5  Filter  bandwidths.  Recall  from  Chapters  II  and  III  that  the  NCDLL  and  MCTL 
BPFs  and  LPFs  have  a  one-sided  noise  bandwidth,  B  Hz,  such  that  B  <  2Wc;  however,  choosing 
too  small  a  value  of  B  increased  the  response  time  of  the  code  tracking  loops  so  as  to  make  simulation 
time  unacceptable.  Taking  these  factors  into  consideration,  we  chose  our  simulation  models  to  have 
filters  with  a  bandwidth  such  that  2Wc/B  =  100  (i.e.,  B  =  0.2  Hz).  Recall  from  chapter  II  that  a 
typical  ratio  for  a  GPS  C/A  code  tracking  loop  is  2Wc/B  «  1000.  The  bandwidths  of  the  LPFs  in 
the  MCEU’s  code  correlator  bank  were  also  chosen  with  bandwidth  B  =  0.2  Hz. 

4.2.6  Code  period  parameter.  Examining  Equation  3  in  Chapter  I,  one  can  see  that  the 
code  correlation  function,  i?c  (^)»  can  be  implemented  by  a  code  multiplier  followed  by  a  BPF  or 
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LPF  with  bandwidth  B  <  2/NTc  Hz,  where  N  is  the  code  period  in  chips  (the  filter  acts  as  an 
integrator).  As  we  have  already  seen,  our  simulation  model  used  a  chip  rate  of  1/Tc  =  10  Hz  with 
filters  of  bandwidth  B  =  0.2  Hz;  this  corresponded  to  a  maximum  code  period  of  iV  =  2”*  -  1  =  63 
chips  (for  m  an  integer),  where  m  =  6  is  the  code  generator’s  shift  register  order  described  in 
Appendix  A.  Choosing  TV  =  63  for  our  simulations  still  allowed  us  to  use  the  large  N  code 
correlation  approximation  of  Equation  5  with  reasonable  accuracy;  however,  the  approximation  is 
not  as  accurate  in  our  case  as  it  is  for  actual  GPS  code  period  parameters  where  N  is  much  greater 
(see  Table  1). 

4.2. 7  Carrier-to-noise  density  ratio  and  loop  SNR.  When  considering  GPS  code  tracking 
loop  performance  in  AWGN,  it  is  important  to  consider  the  received  direct-path  carrier  power  to 
noise  density  ratio  (C/Nq  in  dB-Hz)  and  the  tracking  loop  SNR,  pt  (in  dB).  In  our  case,  the 
received  carrier  power,  G,  equaled  the  transmitted  power,  P ,  since  G  =  Puq  and  we  chose  oq  — 
The  code  tracking  loop’s  SNR  is  determined  by  the  two-sided  PSD  of  the  lowpass  noise  input  to 
the  linear  loop  model  (Figure  8)  and  the  two-sided  loop  noise  bandwidth,  25^  Hz,  where  is 
defined  in  Equation  28. 

For  the  NCDLL,  if  we  assume  TV  »  1  and  an  early-late  offset  of  A  =  1,  the  loop  SNR  is  given 
by  (5) 

For  GPS,  typical  values  of  P/No  range  from  approximately  30  to  45  dB-Hz  while  a  typical  NCDLL 
loop  bandwidth  is  Bx,  =  10  Hz  (4)  (14);  this  gives  a  loop  SNR  between  17  to  32  dB. 

For  the  MCTL,  the  input  noise  PSD  for  the  linear  model  is  defined  over  the  LPF  bandwidth, 
B,  and  can  be  expressed  as  (see  Chapter  III) 


Gne>{f) 


2No{xl  +  xl) 
A2 


(91) 
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where  A  is  the  slope  of  the  MCTL  S-curve  defined  in  Table  2  for  different  signal  parameters 
xo,  xi,  and  a.  Therefore,  the  MCTL  loop  SNR  can  be  written  as 


PL 


P  r>  P 

Gne'2BL  ^  ^^No2Bl 


(for  MCTL) 


(92) 


where  the  loop  SNR  gain,  Gi,  depends  on  the  multipath  signal  parameters  and  is  defined  as 


Gl  =  A^l  {xl  +  X?) 


(93) 


For  all  simulations,  we  chose  P  =  1/2  W  and  the  attenuation  coefficients  as  oo  =  1  and  ai  =  0.5; 
under  our  in-phase  condition  described  in  section  4.2.3,  this  gave  the  multipath  signal  parameters 
as  Xo  =  1  and  xi  =  0.5  (see  Equation  35).  Therefore,  Gi  was  defined  for  different  values  of  a  as 
shown  in  Table  4. 

Table  4.  MRDLL  loop  SNR  gain  (xq  =  1,  =  0,5). 


a 

Gl  (dB) 

0  <  QJ  <  0.5 

12 

a  =  0.5 

8.6 

0.5  <  a  <  1 

2.6 

1  <  a  <  1.5 

2.6 

a  =  1.5 

5 

In  practical  code  tracking  loops,  the  loop  noise  bandwidth  is  usually  much  less  than  the 
BPF  or  LPF  bandwidths  (i.e.,  2Bl  <  B)\  however,  since  our  choice  of  B  was  constrained  to  be 
relatively  small,  satisfying  2Bl  <  B  required  extremely  small  values  of  the  loop  natural  frequency, 
ujn  (since,  from  Equation  28,  Bl  is  proportional  to  a;^).  These  low  values  of  u)n  made  MCTL 
response  slow  and  resulted  in  excessive  simulation  times  and  data  file  sizes;  therefore,  uJn  =  bvr 
rad/sec  was  selected  for  the  MCTL  to  keep  simulations  within  reasonable  limits  while  maintaining 
loop  stability  (the  same  consideration  led  to  choosing  oJn  =  47r  rad/sec  for  the  NCDLL  as  =  Gtt 
rad/sec  made  the  loop  unstable).  However,  this  choice  of  cOn  ~  ^tt  rad/sec  resulted  in  2Bi  ~  20 
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Hz  (obtained  from  Equation  28  for  a  damping  ratio,  ^  =  V^/2),  which  was  100  times  greater  than 
the  MCTL  LPF  bandwidth,  B  =  0.2  Hz;  since  the  input  noise  PSD  for  the  MCTL  linear  model  is 
defined  over  B,  this  LPF  bandwidth  became  the  effective  loop  bandwidth.  Therefore,  the  MCTL 
loop  SNR  for  our  simulations  was 


P«//  =  Gl 


P 

NoB 


(  for  MRDLL  simulation). 


(94) 


For  the  simulations  involving  MRDLL  performance  in  AWGN,  values  of  peff  ranged  from  approx¬ 
imately  20  to  40  dB. 

^.2,8  SimuloMon  parameter  summary.  A  summary  list  of  all  simulation  parameters  is 
presented  below;  any  parameters  not  specifically  covered  in  sections  4.2.2  through  4.2.7  were  chosen 
based  on  designer  preference  and/or  simulation  time  and  data  file  size  constraints: 


•  Sampling  frequency:  fg  =  1000  Hz 

•  Transmitted  signal  power:  P  =  1/2  W 

•  Attenuation  coefficients:  oq  =  1;  cq  =  0.5 

•  Direct-path  propagation  delay:  tq  =  0  sec 

•  Multipath  delay  coefficient:  a  =  O.lfc  (for  k  an  integer  between  0  and  15) 

•  Carrier  phase  (direct-path  and  reflected  signals):  =  9  rad 

•  Local  carrier  phase  error  (MRDLL  baseband  conversion):  (j>e  =0  rad 

•  MRDLL  signal  parameters  (see  Equation  35):  Xo  =  1;  xi  =  0.5  (xi  =  0  for  direct-path  only) 

•  Chip  rate:  1/Tc  =  10  Hz  (i.e.,  1  chip  per  100  samples) 

•  DS/SS  code  period:  JV  =  63  chips/period 

•  Code  Doppler:  negligible 
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•  BPF  bandwidth  (at  IF),  LPF  bandwidth  (two-sided),  and  MCEU  correlator  bank  LPF  band¬ 
width:  B  =  0.2  Hz 

•  NCDLL  and  MCTL  loop  order:  Second-order  with  active  lead-lag  filter  as  described  by 
Equation  24. 

•  NCDLL  natural  frequency:  ujn  —  47r  rad/sec  (from  Equation  28) 

•  MCTL  natural  frequency:  Un  =  Ott  rad/sec  (from  Equation  28) 

•  NCDLL  and  MCTL  damping  ratio:  C  =  V^/2  «  0.707 

•  MCTL  effective  loop  SNR:  20  <  peff  <  40  dB 

•  ALC  gain  threshold  (MRDLL):  KAm&x  =  ^^n/50 

•  Early-late  correlator  spacing  (MRDLL  and  NCDLL):  A  =  1  chip 

•  Delay  spacing  coefficients  (MRDLL):  pk  =  O.U  (for  k  an  integer  between  0  and  15) 


4.3  Data  processing 

All  data  was  processed  in  MATLAB  using  m-files  contained  in  Appendix  B;  for  the  SPW 
simulations,  data  was  first  converted  from  SPW  format  to  MATLAB  mat-file  format  via  SPW’s 
MATLAB  SIGNAL  SINK  block.  Data  processing  focused  mainly  on  two  measurements:  normalized 
root-mean-square  (rms)  code  phase  tracking  error,  and  the  MCEU  estimator  means  and  variances. 


For  each  loop  simulation  run,  the  normalized  rms  tracking  error,  arms  (in  chips),  was  esti¬ 
mated  by  comparing  each  leading  edge  of  the  received  direct-path  DS/SS  code  to  the  corresponding 
leading  edge  of  the  ‘on-time’  replica  code  generated  by  the  code  tracking  loop;  by  representing  the 
number  of  samples  between  the  ith  pair  of  leading  edges  as  di,  we  computed  the  normalized  rms 
tracking  error  estimate  as 


^rms  — 


cfs  M 


chips, 


(95) 
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where  L  was  the  number  of  leading  edge  pairs  present  over  the  given  sample  size.  The  algorithm 
described  by  Equation  95  was  implemented  using  the  leadedge.m  m-file. 

To  measure  MCEU  estimator  performance,  the  MCEU  estimator  outputs  (xg,  xi,  and  3) 
were  sampled  for  a  given  MRDLL  or  MCEU  simulation  run.  The  outputs  were  plotted  versus  time 
using  the  mrdlest.m  m-file;  this  m-file  also  returned  the  sample  means  and  variances  of  xo,xi,  and 
a  using  MATLAB’s  mean  and  cov  functions. 
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44  Simulation  #1:  MCEU  estimator  performance  with  perfect  code  phase  estimation  (MATLAB 
simulation) 

44.1  Objective.  As  seen  in  Chapter  III,  the  MRDLL’s  ability  to  track  out  multipath 
errors  greatly  depends  on  the  ability  of  the  MCEU  to  provide  accurate  signal  parameter  estimates 
to  the  MCTL;  therefore,  investigating  MCEU  estimator  behavior  was  essential  to  evaluating  the  the 
MRDLL  design  for  possible  GPS  multipath  applications.  The  objective  of  this  simulation  was  to 
investigate  the  MCEU  estimator  performance  assuming  zero  MCTL  code  phase  tracking  error  (i.e., 
To  =  To);  by  making  this  assumption,  we  were  able  perform  the  simulation  in  MATLAB  without 
running  the  full  SPW  MRDLL  simulation  model  (which  was  relatively  slow).  The  advantage  of 
using  MATLAB  to  simulate  the  MCEU  estimator  was  that  MATLAB  is  optimized  to  perform  the 
matrix  operations  needed  to  generate  the  estimates  provided  to  the  MCTL. 

4.4. £  Method.  Using  the  mceu.m  m-file,  simulations  were  run  for  different  values  of  a  and 
P/NqB;  each  iteration  of  mceu.m  represented  one  sample  of  the  MCEU  estimator  output.  A  total 
of  100  iterations  were  performed  for  each  value  of  a  and  P/NqB',  then,  for  each  of  these  100-sample 
outputs,  we  estimated  the  sample  mean  and  sample  variance  of  the  ML  estimates  xq,  xi,  and  a. 
In  each  trial,  xq  =  1  and  xi  =  0.5  were  the  true  parameters  except  at  a  =  0,  wherein  xi  =  0  was 
required  (no  multipath  present). 

4.4.3  Results.  In  Figures  26  through  33,  the  mean  values  (taken  over  100  samples)  of  the 
MCEU  ML  estimates  are  plotted  versus  a  for  four  different  values  of  P/NoB;  these  plots  yield  the 
following  observations: 

•  As  expected,  the  mean  values  of  Xo  and  X\  approached  the  actual  signal  parameter  values  as 
P/NoB  increased.  On  the  other  hand,  a  provided  an  accurate  estimate  of  multipath  delay 
only  when  multipath  was  actually  present;  for  the  direct-path  only  case  (i.e.,  a  =  0  and 
Xi  =  0),  a  exhibited  a  significant  bias. 
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P/NoB=:10  dB  (perfect  code  phase  tracking) 


a 


Figure  26.  Simulation  results  showing  the  mean  value  of  MCEU  ML  estimate,  x,  a.t  P/Nq  = 
10  dB,  for  different  values  of  relative  multipath  delay,  a.  (Assumes  perfect  MCTL 
code  phase  tracking). 

•  The  bias  on  3  at  o;  =  0  implied  that,  when  no  multipath  was  present,  the  MCEU  was 
consistently  choosing  the  wrong  Xk  as  its  ML  estimate  and  would  therefore  erroneously  enable 
one  of  the  MCTL’s  multipath  arms.  However,  the  plots  also  show  that  the  bias  on  xq  and  Xi 
at  a  =  0  was  minimal  despite  the  large  bias  on  S;  this  seemingly  contradictory  result  implied 
that  the  effect  on  code  tracking  would  be  minimal,  because  the  signal  in  the  multipath  arm 
would  be  multiplied  by  xi  «  0  (refer  to  Figure  14). 

•  The  ML  estimates  when  the  relative  multipath  delay  was  small  (i.e.,  a  =  0.1)  were  generally 
worse  than  the  ML  estimates  at  greater  relative  multipath  delays. 

•  Correlation  between  the  Xq  and  xi  ML  estimates  is  apparent  (see  Figure  26). 

These  observations  are  explained  theoretically  by  examining  the  statistical  properties  of  the 
16  signal  parameter  estimates,  xj,,  from  which  the  MCEU  chooses  its  ML  estimate,  x.  The  mean 
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P/NoB=15  dB  (perfect  code  phase  tracking) 


a 


Figure  27.  Simulation  results  showing  the  mean  value  of  MCEU  ML  estimate,  x,  at  P/Nq  — 
15  dB,  for  different  values  of  relative  multipath  delay,  a.  (Assumes  perfect  MCTL 
code  phaae  tracking). 


P/NoB=20  dB  (perfect  code  phase  tracking) 


a 


Figure  28.  Simulation  results  showing  the  mean  value  of  MCEU  ML  estimate,  x,  at  P/Nq  = 
20  dB,  for  different  values  of  relative  multipath  delay,  a.  (Assumes  perfect  MCTL 
code  phase  tracking). 
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P/NoB=25  dB  (perfect  code  phase  tracking) 


a 


Figure  29.  Simulation  results  showing  the  mean  value  of  MCEU  ML  estimate,  x,  at  P/Nq  = 
25  dB,  for  different  values  of  relative  multipath  delay,  a.  (Assumes  perfect  MCTL 
code  phase  tracking). 


Figure  30.  Simulation  results  showing  the  mean  value  of  MCEU  ML  estimate,  a,  at  P/iVo  = 
10  dB,  for  different  values  of  relative  multipath  delay,  a.  (Assumes  perfect  MCTL 
code  phase  tracking). 
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P/NoB=15  dB  (perfect  code  phase  tracking) 


a 


Figure  31.  Simulation  results  showing  the  mean  value  of  MCEU  ML  estimate,  S,  at  P/Nq  = 
15  dB,  for  different  values  of  relative  multipath  delay,  a.  (Assumes  perfect  MCTL 
code  phase  tracking). 


P/NoB=20  dB  (perfect  code  phase  tracking) 


a 


Figure  32.  Simulation  results  showing  the  mean  value  of  MCEU  ML  estimate,  a,  at  P/Nq  = 
20  dB,  for  different  values  of  relative  multipath  delay,  a.  (Assumes  perfect  MCTL 
code  phase  tracking). 


P/NoB==25  dB  (perfect  code  phase  tracking) 


a 


Figure  33.  Simulation  results  showing  the  mean  value  of  MCEU  ML  estimate,  3,  at  P/No  = 
25  dB,  for  different  values  of  relative  multipath  delay,  a.  (Assumes  perfect  MCTL 
code  phase  tracking). 

of  xj;  for  ^  0  is  given  by  (see  Chapter  III) 


E{%}  =  E{P  if3k)  Ti}  =  E{P  (/3fc)  [H  (a)  x  +  v]}  =  P  (/34)  H  (a)  x  ,  (96) 


while  the  mean  of  the  k  =  0th  estimate  is 

P(/?o)H(a)x 

0 

When  multipath  is  present,  and  because  P  (^fe)  H  (/J*,)  =  I  for  fc  0,  any  given  estimator  for  fc  7^  0 
is  unbiased  (i.e.,  E  {xj,}  =  x)  when  a  —  the  estimator  for  fc  =  0,  on  the  other  hand,  is  always 
biased  in  the  presence  of  multipath.  For  the  direct-path  only  case,  the  estimator  for  =  0  is 
unbiased  because  P  (;0o )  H  (/Jo)  x  =a:o;  however,  each  of  the  remaining  estimators  for  7^  0  are  also 


i^{xo}  = 
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a  =  0.5;  P/NoB=25dB 
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Figure  34.  Simulation  results  showing  the  mean  values  of  the  MCEU  estimates,  Xjb,  at  P/iVo 
25  dB  when  a  —  0.5.  (Assumes  perfect  MCTL  code  phase  tracking). 


unbiased  since 


P(/3fc)H(0)x  = 


- 

- 

" 

” 

1 

rrti2 

Xq 

Xq 

0 

m22 

0 

0 

(98) 


where  m\2  and  7/122  denote  arbitrary  matrix  elements.  The  fact  that  the  first  column  of  P(/3fe)H(0) 
is  always  [1  O]'^  is  due  to  the  fact  that  'P{/3k)liil3k)  =  I  and  that  the  first  column  of  H(a)  (see 
Equation  70)  is  the  same  regardless  of  the  value  of  a.  The  preceding  discussion  illustrates  an 
important  property  of  the  MCEU  estimator:  when  multipath  is  present,  only  the  Xk  corresponding 
to  /3k  —  a  is  an  unbiased  estimate;  when  no  multipath  is  present,  all  of  the  Xk  are  unbiased 
estimates]  both  cases  are  illustrated  in  Figures  34  through  35. 

To  predict  the  effect  of  the  estimates,  Xfc,  upon  the  MCEU  ML  estimate,  x,  we  calculate  the 
mean  value  of  the  kth  error  vector,  ejb,  as 


E{ek}  =  n{a)x-li{pk)E{xk}  . 


(99) 
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a  ^  0;  P/NoB=25dB 
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Figure  35.  Simulation  results  showing  the  mean  values  of  the  MCEU  estimates,  x*,,  at  P/Nq  — 
25  dB  when  a  =  0  and  xi  =  0  (i.e.,  no  multipath).  (Assumes  perfect  MCTL  code 
phase  tracking). 


From  the  preceding  paragraph  we  see  that,  when  no  multipath  is  present,  E  {e*}  =  0  when  /J*  =  a. 
Therefore,  since  Ek  =  efC-^efc  =  0  maximizes  the  likelihood  function  in  Equation  89,  the  MCEU 
will  tend  to  choose  the  estimates  corresponding  to  /3fc  =  a  as  its  ML  estimates  (note,  however,  that 
E{Ek}  7^  0).  For  the  direct-path  only  case,  E  {cfc}  =  0  for  all  values  of  k,  which  makes  any  choice 
of  pk  acceptable  as  far  as  the  MCEU  is  concerned;  this  explains  why  the  a  estimate  can  exhibit 
a  large  bias  while  the  xq  and  Xi  estimates  remain  accurate  (as  demonstrated  in  the  simulation 
results).  Therefore,  the  MRDLL  can  function  in  both  the  direct-path  and  multipath  environments. 

Now,  we  will  examine  the  covariance  matrix,  W,  for  x*,,  which  can  be  shown  to  be  (6) 


W  =  a3[HT  (/30C-iH(^k)] 


(100) 


and  is  2  X  2  symmetric.  The  Xj,  estimator  variance  (taken  from  the  diagonal  terms  of  W)  is  plotted 


for  different  values  of  P/NoB  in  Figures  36  through  43  along  with  the  sample  variances  of  the 
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P/NoB=10  dB  (perfect  code  phase  tracking) 


Figure  36.  Simulation  results  showing  the  variance  of  MCEU  ML  estimate,  ®o,  at  P/Nq  -  10  dB, 
for  different  values  of  relative  multipath  delay,  a.  (Assumes  perfect  MCTL  code  phase 
tracking). 


simulated  ML  estimates.  Note  that  the  variance  of  the  estimates  at  a  =  0.1  is  significantly  greater 
than  that  for  the  other  estimates;  this  explains  why  the  estimates  at  small  relative  multipath  delay 
(a  —  0.1)  were  generaUy  not  as  accurate  as  those  for  other  values  of  delay.  In  addition,  calculation 
of  W  for  different  values  of  fik  reveals  the  presence  of  non-zero  diagonal  terms  for  0  <  /3fc  <  1, 
which  implies  a  cross-correlation  exists  between  x^k  and  x\k\  this  explains  the  correlation  observed 
in  Figure  26. 
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0.35 


P/NoB=10  dB  (perfect  code  phase  tracking) 


a 


Figure  37.  Simulation  results  showing  the  variance  of  MCEU  ML  estimate,  xi,  at  P/Nq  =  10  dB, 
for  different  values  of  relative  multipath  delay,  a.  (Assumes  perfect  MCTL  code  phase 
tracking). 


P/NoB=15  dB  (perfect  code  phase  tracking) 


Figure  38.  Simulation  results  showing  the  variance  of  MCEU  ML  estimate,  xo,  at  P/Nq  =  15  dB, 
for  different  values  of  relative  multipath  delay,  a.  (Assumes  perfect  MCTL  code  phase 
tracking). 
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Figure  39,  Simulation  results  showing  the  variance  of  MCEU  ML  estimate,  xi,  at  P/Nq  =  15  dB, 
for  different  values  of  relative  multipath  delay,  a,  (Assumes  perfect  MCTL  code  phase 
tracking). 


P/NoB=20  dB  (perfect  code  phase  tracking) 


Figure  40,  Simulation  results  showing  the  variance  of  MCEU  ML  estimate,  xo^  P/Nq  —  20  dB, 
for  different  values  of  relative  multipath  delay,  a,  (Assumes  perfect  MCTL  code  phase 


P/NoB-20  dB  (perfect  code  phase  tracking) 


Figure  41.  Simulation  results  showing  the  variance  of  MCEU  ML  estimate,  xi,  at  P/Nq  =  20  dB, 
for  different  values  of  relative  multipath  delay,  a.  (Assumes  perfect  MCTL  code  phase 
tracking). 


Figure  42.  Simulation  results  showing  the  variance  of  MCEU  ML  estimate,  xq,  at  P/iVo  —  25  dB, 
for  different  values  of  relative  multipath  delay,  a.  (Assumes  perfect  MCTL  code  phase 
tracking). 
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4.5  Simulation  #2:  NCDLL  vs  MRDLL  steady-state  tracking  performance  in  the  presence  of 
multipath  without  AWGN 

4.5.1  Objectives.  The  effect  of  multipath  interference  on  the  NCDLL  (neglecting  noise) 
was  discussed  in  Chapter  II;  we  saw  that  the  reflected  signal  introduced  a  steady-state  code  phase 
tracking  error  in  the  NCDLL.  Then,  in  Chapter  III  we  claimed  that  the  MCTL  could  remove 
this  tracking  error  by  effectively  removing  the  reflected  signal’s  contribution  to  the  loop’s  S-curve. 
However,  the  analysis  in  Chapter  III  assumed  that  the  MCTL  received  perfect  signal  estimates  from 
the  MCEU;  furthermore,  the  MCEU  simulation  in  section  4.4  assumed  zero  code  phase  tracking 
error.  Therefore,  in  order  to  view  MRDLL  behavior  under  conditions  where  these  ideal  conditions 
no  longer  apply,  computer  simulation  was  necessary;  the  objectives  of  this  simulation  were  to 

1.  Investigate  the  effects  of  multipath  interference  on  the  MCTL  steady-state  tracking  error  with 
imperfect  MCEU  signal  parameter  estimation  (neglecting  noise  effects). 

2.  Investigate  the  effects  of  multipath  interference  on  the  MCEU  estimator  when  code  phase 
timing  errors  are  present  (without  AWGN). 

3.  Demonstrate  that  the  MRDLL  exhibits  reduced  steady-state  code  phase  tracking  error  when 
compared  to  the  NCDLL  in  the  presence  of  multipath  without  AWGN. 

4.5.2  Method.  Simulations  were  run  for  both  loops  at  different  values  of  a  >  0  (the 
direct-path  only  simulation  at  a  =  0  is  presented  in  the  next  section).  Since  we  were  interested 
in  loop  performance  after  acquisition,  both  loops  were  allowed  to  reach  zero  steady-state  tracking 
error  before  introduction  of  the  multipath  signal;  in  the  MRDLL,  this  was  accomplished  by  feeding 
the  MCTL  perfect  estimates  at  the  start  of  each  run,  then  switching  to  the  MCEU  estimates. 
After  the  introduction  of  the  multipath  signal,  each  loop  was  allowed  to  return  to  steady-state;  the 
normalized  rms  steady-state  timing  error  was  then  estimated  over  the  last  100,000  output  samples 
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using  the  method  described  in  section  4.3.  The  sample  mean  and  sample  variance  of  the  MCEU 
estimates  at  steady-state  were  also  calculated  over  the  last  100,000  output  samples. 

4.5.3  Results.  The  normalized  rms  steady-state  tracking  error  resulting  from  the  NCDLL 
and  MRDLL  simulations  is  plotted  versus  a  in  Figure  44.  This  plot  shows  that: 

•  The  MRDLL  exhibits  a  smaller  tracking  error  than  the  NCDLL  for  all  values  of  a;  also  note 
that  the  improvement  is  most  significant  at  values  of  a  which  introduce  the  greatest  error  in 
the  NCDLL. 

•  Since  we  no  longer  have  perfect  MCEU  estimation,  the  MRDLL  exhibits  a  non-zero  steady- 
state  error. 

•  The  MRDLL  steady-state  error  peaks  at  a  =  0.6,  much  in  the  same  way  the  NCDLL  error 
peaks  between  a  =  0.7  and  a  =  0.8. 

•  The  NCDLL  simulation  results  agree  almost  exactly  with  the  multipath  tracking  error  pre¬ 
dicted  by  Equation  32. 

MCEU  estimator  performance  is  plotted  in  Figures  45  through  49;  in  addition,  several  plots 
showing  the  MCEU  output  over  time  for  selected  values  of  a  are  presented  in  Figures  50  through 
55.  Plots  showing  the  MCTL  initial  transient  response  after  introduction  of  the  multipath  signal 
are  also  shown  in  Figures  56  through  63.  (Note  that  the  time  scaling  on  the  MCEU  plots  is  based 
on  the  simulation  parameter  of  Tc  =  0.1  seconds;  keep  in  mind  that  actual  GPS  time  would  scale 
our  time  axis  by  a  factor  on  the  order  of  10“^  seconds  since,  assuming  C/A  code,  the  actual  GPS 
code  chip  period  is  approx  10"®).  From  the  simulation  results,  observe  that: 

•  As  expected,  the  imperfect  code  phase  estimation  introduced  a  bias  into  the  MCEU  ML  signal 
estimates,  xq  and  Xi. 

•  For  the  xi  estimate,  the  bias  was  relatively  large  at  a  =  0.6,  corresponding  to  the  peak 
steady-state  timing  error  observed  in  Figure  44. 
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MRDLLvs  NCDLL 


a 


Figure  44.  Simulation  results  comparing  MRDLL  and  NCDLL  rms  steady-state  tracking  error 
performance  for  different  values  of  relative  multipath  delay,  a. 

•  The  variance  of  the  xq  and  xi  ML  estimates  tended  to  increase  with  increasing  a;  this  is  the 
opposite  of  what  was  observed  in  Simulation  #1. 

•  There  were  significant  increases  at  a  =  0.6  in  the  variances  of  xo  and  3;  in  fact  the  3  ML 
estimate  when  a  =  0.6  was  the  only  3  with  a  non-zero  variance  at  steady-state. 

•  The  MCEU  outputs  display  an  oscillatory  behavior  at  steady-state,  which  is  likely  due  to  the 
presence  of  adaptation  in  the  MCTL  loop  filter  (via  the  ALC). 

The  most  significant  result  from  these  simulations  is  that,  despite  a  non-zero  steady-state 
tracking  error,  the  MRDLL  demonstrates  a  dramatic  improvement  over  the  NCDLL  in  code  phase 
tracking  error.  The  MRDLL  non-zero  steady  state  error  exists  because  of  the  interdependence 
between  the  MCEU  estimates  and  the  MCTL  code  phase  estimate;  the  code  phase  error,  affects 
the  accuracy  of  the  x  estimate  and  vice-versa,  which  is  evident  in  the  oscillatory  nature  of  the 
MCEU  output.  Furthermore,  the  peaks  in  the  steady-state  error  at  a  =  0.6  may  be  due  to  the  fact 
that  the  magnitude  of  the  peak  transient  error  varied  with  a,  and  was  greatest  at  a  =  0.6;  perhaps 
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MRDLL  multipath  only  simulation 


a 


Figure  45.  Simulation  results  showing  the  mean  value  of  MCEU  ML  estimate,  x  for  different 
values  of  relative  multipath  delay,  Q;  (neglects  noise  effects). 


MRDLL  multipath  only  simulation 
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Figure  46.  Simulation  results  showing  the  mean  value  of  MCEU  ML  estimate,  a,  for  different 
values  of  relative  multipath  delay,  a  (neglects  noise  effects). 
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MRDLL  multipath  only  simulation 


Figure  47.  Simulation  results  showing  the  variance  of  MCEU  ML  estimate,  xq,  for  different  values 
of  relative  multipath  delay,  a  (neglects  noise  effects). 


Figure  48.  Simulation  results  showing  the  variance  of  MCEU  ML  estimate,  xi,  for  different  values 
of  relative  multipath  delay,  a  (neglects  noise  effects). 
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Figure  49. 


Figure  50. 
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Simulation  results  showing  the  variance  of  MCEU  ML  estimate,  a,  for  different  values 
of  relative  multipath  delay,  a  (neglects  noise  effects). 


alpha=0.3 


Simulation  results  showing  MCEU  estimator  output,  x,  over  time  for  a  =  0.3  (no 
noise  L 
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Figure  55.  Simulation  results  showing  MCEU  estimator  output,  a,  over  time  for  a  =  0.7  (no 
noise). 


in  this  case,  an  error  threshold  was  exceeded  that  prevented  the  code  phase  and  signal  estimates 
from  attaining  more  accurate  steady  state  values  (compare  the  MCEU  outputs  at  steady  state  for 
a  =  0.6  to  those  for  a  =  0.3  and  a  =  0.7). 
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Figure  56.  Simulation  results  showing  MRDLL  transient  code  phase  timing  error  for  a  =  0.1  and 
a  —  0.2  (no  noise). 
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Figure  57.  Simulation  results  showing  MRDLL  transient  code  phase  timing  error  for  a  =  0.3  and 
a  =  0.4  (no  noise). 
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Figure  58.  Simulation  results  showing  MRDLL  transient  code  phase  timing  error  for  a  —  0.5  and 
a  =  0.6  (no  noise). 
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Figure  59,  Simulation  results  showing  MRDLL  transient  code  phase  timing  error  for  a  =  0.7  and 
a  =  0,8  (no  noise). 
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Figure  60.  Simulation  results  showing  MRDLL  transient  code  phase  timing  error  for  a  -  0.9  and 
a  =  1.0  (no  noise). 
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Figure  61.  Simulation  results  showing  MRDLL  transient  code  phase  timing  error  for  a  -  1.1  and 
CL  —  1.2  (no  noise). 
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Figure  62.  Simulation  results  showing  MRDLL  transient  code  phase  timing  error  for  a  —  1.3  and 
a  —  1.4  (no  noise). 
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Figure  63.  Simulation  results  showing  MRDLL  transient  code  phase  timing  error  for  a  =  1.5  (no 


J^.G  Simulation  #3:  MRDLL  vs  NCDLL  steady-state  tracking  performance  in  a  direct-path  only 

environment  without  AWGN, 

4.6.1  Objective.  As  we  saw  from  the  previous  simulation,  the  MRDLL  exhibits  a  signifi¬ 
cant  improvement  in  multipath  code-phase  tracking  error  over  the  NCDLL;  however,  the  MRDLL 
must  also  be  able  to  track  the  GPS  code  phase  when  no  reflected  signal  is  present  to  have  any  prac¬ 
tical  application.  The  objective  of  this  simulation  was  to  examine  MRDLL  steady-state  tracking 
error  when  no  multipath  is  present,  and  compare  it  to  that  of  the  NCDLL. 

4.6.2  Method.  This  simulation  was  identical  to  Simulation  #2;  the  only  difference  was 
that  a  direct-path  signal  (i.e.,  a  =  0  and  xi  =  0)  was  introduced  after  each  loop  entered  the  initial 
steady-state.  The  rms  tracking  error  and  MCEU  estimator  means  and  variances  were  calculated 
over  100,000  samples  as  in  Simulation  #2. 

4.6.3  Results.  Simulation  results  are  presented  in  Table  5.  Note  from  the  table  that 
the  MRDLL  exhibits  a  small  but  non-zero  steady-state  code  tracking  error  while  for  the  NCDLL, 
aT,rn3  =  0;  therefore,  the  NCDLL  is  a  better  choice  for  code  phase  tracking  in  the  direct-path  only 
environment.  The  reason  for  this  is  simple:  the  MRDLL  must  estimate  the  amplitude  (or  existence) 
of  the  reflected  signal,  whereas  the  NCDLL  design  already  assumes  that  no  multipath  is  present. 

Observing  the  MCEU  ML  estimates,  note  that  despite  the  introduction  of  code  phase  tracking 
error,  xq  «  1  and  Xi  «  0;  however,  the  a  estimate  shows  a  significant  bias  and  a  relatively  high 
variance.  These  results  are  similar  to  the  direct-path  results  of  the  MATLAB  MCEU  simulation  in 
section  4.4. 
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Table  5. 


MRDLL  vs  NCDLL  in  direct-path  only  environment  (xo  =  1,  xi  =  0,  a 
results. 


Parameter 

MRDLL 

NCDLL 

arms  (chips) 

0.0036 

0.0000 

Xo  (mean) 

0.9871 

N/A 

xi  (mean) 

-0.0158 

N/A 

a  (mean) 

1.2376 

N/A 

Xo  (variance) 

2.6432  X  10“** 

N/A 

Xi  (variance) 

.0056 

N/A 

a  (variance) 

.0634 

N/A 

0):  simulation 
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^.7  Simulation  #4:  MRDLL  steady-state  tracking  performance  in  the  presence  of  multipath  with 
AWGN 

4.7.1  Objectives.  Now  that  we  have  demonstrated  the  MRDLL ’s  ability  to  significantly 
reduce  code  phase  tracking  errors  due  to  multipath,  we  desire  to  investigate  the  effects  of  AWGN 
on  MRDLL  tracking  performance  in  the  multipath  environment.  Recall  that  the  noise  analysis 
and  models  in  Chapter  III  assumed  perfect  MCEU  estimation  and  MCTL  code  phase  tracking  in 
the  S-curve’s  linear  operating  region.  In  practical  tracking  loops,  linear  tracking  occurs  only  when 
the  SNR  exceeds  a  certain  threshold;  computer  simulation  is  a  useful  tool  for  both  demonstrating 
this  threshold  effect  as  well  as  confirming  theoretical  predictions  based  on  the  assumption  of  linear 
tracking  and  perfect  estimation.  The  objectives  of  this  simulation  were  to: 

1.  Investigate  MRDLL  non-linear  behavior  in  the  presence  of  multipath  with  AWGN,  demon¬ 
strating  the  SNR  threshold  effect. 

2.  Compare  MRDLL  simulated  linear  tracking  performance  in  AWGN  without  perfect  MCEU 
estimation  to  that  predicted  by  Equation  61  for  perfect  estimation. 

3.  Evaluate  the  effects  of  non-linear  tracking  on  the  MCEU  estimator. 

4.7.2  Method.  Simulations  were  run  for  different  values  of  MRDLL  loop  SNR,  peff  (in 
dB)  and  five  different  values  of  a  =  (0.4,  0.5,  0.7,  1.1,  1.5);  since  we  were  investigating  linear  versus 
non-linear  behavior  and  since  the  MCTL  linear  region  varied  with  a,  we  chose  one  value  of  a  from 
each  of  the  five  linear  regions  described  in  Table  2.  As  in  the  previous  loop  simulations,  the  MRDLL 
was  allowed  to  reach  steady-state  prior  to  introduction  of  the  multipath  signal;  once  the  MCTL 
returned  to  steady-state,  arms  and  MCEU  estimator  mean  and  variance  were  calculated  over  the 
last  20,000  output  samples. 

4.7.3  Results.  Simulation  results  plotting  arms  versus  pe/f  are  shown  in  Figures  64 
through  68;  this  data  yields  the  following  observations: 
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a  =  0.4 


Figure  64.  Simulation  results  showing  estimated  steady-state  tracking  error,  cTrms  (in  chips),  for 
different  values  of  effective  MCTL  loop  SNR,  Peff  (in  dB)  when  a  =  0.4. 

•  For  each  different  value  of  a,  there  there  is  an  obvious  loop  SNR  threshold  (in  dB)  at  which 
the  loop  tends  toward  linear  tracking  behavior  (denoted  by  the  dashed  line). 

•  The  SNR  threshold  is  not  constant,  but  depends  on  the  value  of  a. 

•  Simulation  results  in  the  linear  region  closely  agree  with  those  predicted  by  Equation  61. 

Figures  69  through  78  plot  the  means  of  the  MCEU  ML  estimates,  along  with  the  results 
from  Simulation  #1  for  comparison;  (the  variances  of  the  ML  estimates  are  included  in  Appendix 
C).  Also,  Figures  79  through  82  represent  typical  MCEU  outputs  versus  time  for  both  low  and 
high  values  of  />«//•  From  the  MCEU  estimator  plots,  observe  that: 

•  The  ML  estimates  exhibit  a  significant  bias  in  the  non-linear  tracking  region,  and  the  point 
at  which  linear  tracking  begins  is  clearly  evident  from  the  plots  for  xo  and  xi;  for  a  given  a, 
note  how  this  point  corresponds  to  the  SNR  threshold  in  Figures  64  through  68. 
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Figure  65.  Simulation  results  showing  estimated  steady-state  tracking  error,  (Jrma  (hi  chips),  for 
dilferent  values  of  effective  MCTL  loop  SNR,  Pe//  (in  <iB)  when  a  =  0.5. 


Figure  66.  Simulation  results  showing  estimated  steady-state  tracking  error,  (Trma  (in  chips),  for 
different  values  of  effective  MCTL  loop  SNR,  peff  (in  dB)  when  a  =  0.7. 


Loop  SNR,  peff  (dB) 

Figure  67.  Simulation  results  showing  estimated  steady-state  tracking  error,  arms  (in  chips),  for 
different  values  of  effective  MCTL  loop  SNR,  peff  (in  dB)  when  a  =  1.1. 


a=  1.5 


Loop  SNR,  peH  (dB) 

Figure  68.  Simulation  results  showing  estimated  steady-state  tracking  error,  arms  (in  chips),  for 
different  values  of  effective  MCTL  loop  SNR,  petf  (in  dB)  when  a  =  1.5. 


•  Above  the  SNR  threshold,  the  code  phase  tracking  error  maintained  linear  performance  even 
when  signal  estimates  were  relatively  poor. 

•  As  expected,  code  phase  estimation  errors  generally  introduce  an  additional  bias  into  the 
mean  ML  estimates  (in  addition  to  that  introduced  by  AWGN);  this  is  seen  by  comparing  the 
data  points  corresponding  to  perfect  code  phase  estimation  with  AWGN  (Simulation  #1). 

Simulation  results  show  that  code  phase  tracking  in  the  linear  regions  followed  theoretical 
predictions  very  closely;  the  non-linear  behavior  of  the  MRDLL,  however,  is  difficult  to  predict  by 
theory,  but  one  can  hypothesize  that  the  fact  that  the  loop  SNR  threshold  varied  with  a  follows 
from  the  fact  that  the  MCTL  linear  operating  region  varies  with  a  as  well  (although  the  complete 
tracking  region  of  positive  slope,  termed  the  affine  region,  does  not  depend  upon  a  and  is  always 
|5|  <  0.5).  For  instance,  note  that  the  threshold  for  a  =  0.4  entered  the  linear  tracking  region  at  a 
loop  SNR  approximately  5  dB  higher  than  that  for  a  =  0.5;  this  corresponds  to  the  fact  that  the 
MCTL’s  linear  tracking  region  when  a  =  0.4  (in  terms  of  normalized  tracking  error,  6)  is  one-fifth 
that  of  a  =  0.5  (see  Table  2).  Also  note  that  the  simulations  at  a  =  1.1  and  a  =  1.5  generally 
demonstrated  a  lower  SNR  threshold  than  the  other  three  values  of  cc;  this  is  most  likely  due  to  the 
fact  that  the  MCTL  at  these  values  of  a  has  a  relatively  large  linear  operating  region  (6  =  ±0.4 
chips  and  8  =  ±0.5  chips,respectively)  and  the  fact  that  transient  and  steady-state  tracking  errors 
are  lower  (as  seen  in  Simulation  #2).  Finally,  the  fact  that  linear  phase  tracking  performance 
(above  the  SNR  threshold)  was  achieved  even  when  signal  estimates  were  poor,  shows  that  SNR 
and  not  estimator  performance  is  the  limiting  factor  for  linear  loop  performance.  This  particular 
result  demonstrates  that  this  simulation  would  be  an  excellent  simulation  to  perform  with  real  GPS 
signal  parameters,  as  it  would  help  determine  limits  on  the  MRDLL  operational  environment  in 
terms  of  received  signal  SNR. 
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Figure  69.  Simulation  results  showing  MCEU  ML  estimate,  x  for  different  values  of  effective 
MCTL  loop  SNR,  pe/f  (in  dB)  when  a  =  0.4. 
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Figure  70.  Simulation  results  showing  MCEU  ML  estimate,  a  for  different  values  of  effective 
MCTL  loop  SNR,  peff  (in  dB)  when  a  =  0.4. 
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Figure  71.  Simulation  results  showing  MCEU  ML  estimate,  x  for  different  values  of  effective 
MCTL  loop  SNR,  peff  (in  dB)  when  a  =  0.5. 
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Figure  72.  Simulation  results  showing  MCEU  ML  estimate,  a  for  different  values  of  effective 
MCTL  loop  SNR,  Peff  (in  dB)  when  a  =  0.5. 
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Figure  73,  Simulation  results  showing  MCEU  ML  estimate,  x  for  different  values  of  effective 
MCTL  loop  SNR,  pejf  (in  dB)  when  a  =  0.7. 
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Figure  74.  Simulation  results  showing  MCEU  ML  estimate,  a  for  different  values  of  effective 
MCTL  loop  SNR,  peff  (in  dB)  when  a  —  0.7. 
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Figure  75.  Simulation  results  showing  MCEU  ML  estimate,  x  for  different  values  of  effective 
MCTL  loop  SNR,  />«//  (in  dB)  when  a  =  1.1. 
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Figure  76.  Simulation  results  showing  MCEU  ML  estimate,  a  for  different  values  of  effective 
MCTL  loop  SNR,  Peff  (in  dB)  when  a  =  1.1. 
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Figure  77.  Simulation  results  showing  MCEU  ML  estimate,  x  for  different  values  of  effective 
MCTL  loop  SNR,  pe/f  (in  dB)  when  a  =  1.5. 


MRDLL  SNR  simulation,  a  =1.5 
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Figure  78.  Simulation  results  showing  MCEU  ML  estimate,  a  for  different  values  of  effective 
MCTL  loop  SNR,  peff  (in  dB)  when  a  =  1.5. 
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V.  Conclusions  and  Recommendations 


5. 1  Overview 

This  thesis  proposed  a  new  code  tracking  loop,  the  modified  RAKE  delay-lock  loop  (MRDLL), 
for  GPS  multipath  mitigation.  Prior  to  introducing  the  proposed  loop,  the  GPS  multipath  problem 
was  examined  by  observing  the  effects  of  multipath  interference  on  the  the  standard  noncoherent 
delay-lock  loop  (NCDLL),  which  is  typically  used  by  GPS  receivers.  Once  the  effects  of  the  multi- 
path  problem  were  understood,  we  proposed  a  solution  in  the  form  of  the  MRDLL;  this  was  done  by 
performing  a  theoretical  analysis  of  the  MRDLL’s  three  main  components:  the  multiple-correlator 
tracking  loop  (MCTL),  the  adaptive  loop  controller  (ALC),  and  the  multiple-correlator  estima¬ 
tion  unit  (MCEU).  Finally,  after  theoretical  predictions  indicated  the  MRDLL’s  ability  to  mitigate 
multipath  effects  on  code  tracking,  computer  simulations  were  performed  to  verify  theoretical  pre¬ 
dictions  and  investigate  loop  behavior  in  the  absence  of  simplifying  assumptions.  This  chapter 
provides  a  summary  of  the  research  results,  and  includes  recommendations  for  further  research. 

5.2  NCDLL  results 

We  began  the  NCDLL  analysis  by  examining  the  equations  that  govern  the  loop’s  behavior; 
in  doing  so,  we  demonstrated  the  ability  of  the  NCDLL  to  track  the  code  phase  of  the  incoming 
GPS  spreading  code  signal  in  a  direct-path  only  environment.  The  analysis  showed  that  under 
direct-path  conditions,  the  NCDLL  tended  to  operate  at  the  point  on  its  discriminator  output 
curve,  or  ‘S-curve’,  where  the  steady-state  code  phase  error  was  zero.  The  NCDLL  linear  model 
was  also  developed;  this  model  became  a  valuable  tool  in  the  subsequent  MRDLL  analysis  and 
simulation. 

Next,  the  effect  of  multipath  interference  on  GPS  code  tracking  via  its  influence  on  the  NCDLL 
steady-state  tracking  point  was  investigated;  the  significance  of  code  tracking  error  was  emphasized 
by  pointing  out  that  a  seemingly  small  error  in  code  phase  estimation  can  translate  into  substantial 
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GPS  ranging  error.  The  NCDLL  multipath  analysis,  performed  by  neglecting  noise  effects  and 
focusing  on  the  effects  of  a  single  reflected  signal,  showed  that  the  reflected  signal  introduced  a 
bias  in  the  NCDLL’s  steady-state  tracking  point;  this  was  caused  by  the  introduction  of  additional 
cross-terms  and  delay  shifts,  which  formed  a  composite  S-curve  no  longer  centered  about  the  zero- 
phase  tracking  point.  In  addition,  this  composite  S-curve  was  demonstrated  to  be  a  function  of 
the  received  multipath  signal  parameters  (i.e.,  amplitudes  of  and  relative  delay  between  direct-path 
and  reflection);  this  implied  that,  for  a  GPS  multipath  mitigation  scheme  to  be  effective,  accurate 
gain  and  pha.se  estimates  would  be  required  to  remove  the  contribution  of  the  reflected  signal.  The 
analysis  also  revealed  that  multipath  errors  were  greatest  when  the  reflected  signal  was  in  carrier 
phase  with  the  direct-path  signal. 

5.3  MRDLL  results 

5.3.1  Theoretical  analysis.  After  demonstrating  the  NCDLL  code  phase  tracking  degra¬ 
dation  caused  by  multipath  interference,  we  proposed  a  solution  to  the  problem  in  the  form  of 
the  MRDLL;  the  subsequent  analysis  showed  that  the  MRDLL  had  the  capability  to  remove  the 
contribution  of  a  reflected  signal  and  track  the  code  phase  with  zero  steady-state  error  in  the  pres¬ 
ence  of  multipath.  The  MRDLL  was  a  modified  version  of  a  code-division  multiple  access  (CDMA) 
tracking  loop  (the  RAKE  delay-lock  loop)  introduced  in  (7);  the  main  difference  was  that  we  de¬ 
signed  the  MRDLL  to  mitigate  the  effects  of  small  relative  multipath  delays  (less  than  1.5  code 
chips),  whereas  the  ROLL  was  designed  for  delays  corresponding  to  integer  multiples  of  the  code 
chip  period. 

The  MRDLL’s  tracking  loop  component,  the  MCTL,  is  a  design  similar  to  that  of  the  NCDLL; 
both  loops  rely  on  early-late  code  correlation  to  track  the  received  signal’s  code  phase  and  provide 
the  receiver  with  an  ‘on-time’  replica  code.  However,  the  MCTL  does  not  perform  a  squaring 
operation  on  the  signal  since  phase  information  is  needed  to  perform  signal  estimation;  this  signal 
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estimation  was  shown  to  be  crucial  to  the  MRDLL’s  ability  to  determine  the  contribution  of  the 
reflected  signal  and  remove  its  effects.  Another  key  difference  between  the  NCDLL  and  MRDLL  is 
that  the  MRDLL  consists  of  multiple  tracking  arms;  one  corresponding  to  the  direct-path  signal 
component  and  the  others  arms  to  candidate  delays  of  a  reflected  signal. 

Assuming  perfect  signal  parameter  estimates  from  the  MCEU,  analysis  showed  that  the 
MCTL’s  linear  tracking  point  was  formed  (as  in  the  NCDLL  multipath  case)  from  a  summation  of 
delay-shifted  and  attenuated  operating  curves;  however,  unlike  the  NCDLL,  the  operating  curves 
were  combined  (using  signal  estimates  from  the  MCEU),  to  yield  a  composite  S-curve  possessing  a 
steady-state  tracking  point  of  zero.  Therefore,  the  MCTL  (given  perfect  signal  estimation),  is  able 
to  successfully  track  the  GPS  code  phase  in  a  multipath  environment. 

Analysis  of  the  MCTL’s  S-curve  and  the  resulting  linear  model  demonstrated  a  very  interest¬ 
ing  property  of  the  MCTL:  the  MRDLL’s  linear  operating  region  varies  with  the  multipath  signal 
signal  parameters.  This  realization  led  to  the  development  of  the  ALC,  which  was  introduced  to 
allow  the  loop  designer  to  fix  the  linearized  dynamic  parameter  specifications  of  the  loop  such  as 
natural  frequency,  damping  ratio,  and  noise  bandwidth.  The  ALC  was  designed  as  a  loop  filter 
preceded  by  an  adaptive  gain  which  adjusted  its  value  according  to  estimates  received  from  the 
MCEU,  once  again  demonstrating  the  importance  of  signal  estimation  to  optimal  tracking  loop 
performance.  The  ALC  design  was  not  present  in  the  RDLL  and  was  unique  to  the  MRDLL 
design. 

The  MCEU  analysis  was  performed  by  assuming  perfect  code  phase  estimates  at  the  MCEU’s 
16  correlator  bank  inputs;  these  correlators,  in  effect,  sample  the  code  correlation  function  at  16  dif¬ 
ferent  values  of  relative  multipath  delay.  The  MCEU  estimator  uses  these  correlator  measurements 
and  knowledge  of  the  noise  output  correlation  matrix  to  generate  16  candidate  maximum-likelihood 
(ML)  signal  estimates,  each  determined  assuming  a  particular  value  of  relative  multipath  delay; 
then,  the  MCEU  chooses  the  most-likely  of  these  16  estimates  via  examination  of  the  resulting 
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prediction  error  vectors.  As  already  discussed,  these  estimates  are  used  by  the  MCTL  to  remove 
the  tracking  error  introduced  by  the  reflected  signal. 

5,3.2  Computer  simulations.  Computer  simulations  were  performed  to  verify  the  results 
predicted  by  theoretical  analysis,  as  well  as  to  examine  MRDLL  behavior  in  the  absence  of  sim¬ 
plifying  assumptions  such  as  perfect  signal  estimation  and  linear  operation.  The  initial  intent  of 
the  simulations  was  to  simulate  (as  close  as  possible)  real-world  GPS  conditions;  however,  the  Sig¬ 
nal  Processing  Workstation  (SPW)  software  simulation  package  proved  to  be  prohibitively  slow  in 
processing  closed  loop  simulations  of  the  MRDLL.  Therefore,  the  focus  of  the  simulations  became 
to  demonstrate  general  MRDLL  behavior  and  tendencies  in  a  variety  of  multipath  and  direct-path 
environments;  the  results  provided  valuable  insight  into  the  application  of  the  MRDLL  to  the  GPS 
multipath  problem.  It  should  be  mentioned  that  the  MATLAB  software  performed  extremely  well 
(in  terms  of  simulation  time)  in  the  MCEU  simulation  assuming  perfect  MCTL  code  phase  esti¬ 
mates;  this  was  because  MATLAB  is  optimized  for  matrix  operations  that  are  used  by  the  MCEU 
estimator  (SPW’s  matrix  algorithms  are  not  even  remotely  as  efficient). 

The  first  simulation  (using  MATLAB)  investigated  MCEU  performance  in  AWGN  assuming 
perfect  code  phase  estimation.  The  results  demonstrated  the  ability  of  the  MCEU  to  accurately 
estimate  the  received  signal  parameters  in  both  direct-path  and  multipath  environments.  It  was 
seen  that,  when  multipath  was  present,  the  MCEU  estimate  corresponding  to  the  received  relative 
multipath  delay  was  the  only  unbiased  estimate  out  of  the  16  possible;  this  meant  the  MCEU  tended 
to  send  this  estimate  as  its  ML  estimate  to  the  MCTL.  For  the  direct-path  only  environment, 
however,  it  was  observed  that  the  MCEU  tended  to  choose  the  estimate  that  did  not  correspond  to 
zero  relative  delay,  but  still  provided  accurate  signal  amplitude  estimates;  this  result  was  due  to  the 
fact  that,  when  no  reflection  was  present,  all  of  the  16  possible  estimates  were  unbicused  estimators. 
The  simulation  also  revealed  that  the  estimator  performance  was  worse  when  the  reflected  signal 


99 


exhibited  a  small  relative  multipath  delay;  this  was  because  the  estimator  variance  depended  on 
the  delay  value  and  was  greatest  when  the  relative  delay  was  small. 

The  second  and  third  simulations  used  the  SPW  environment  to  compare  MRDLL  vs  NCDLL 
steady-state  code  phase  tracking  performance  in  the  presence  of  multipath  interference.  The  results 
showed  that,  as  predicted  by  the  theoretical  analysis,  the  MRDLL  exhibited  a  significant  improve¬ 
ment  in  steady-state  tracking  error  over  the  NCDLL  when  multipath  was  present]  the  presence  of  a 
non-zero  steady-state  error,  however,  indicated  the  effects  of  imperfect  estimation  on  the  MRDLL. 
In  a  direct-path  only  environment,  it  was  the  NCDLL  that  demonstrated  the  best  tracking  per¬ 
formance  as  it  had  zero  steady-state  error  compared  to  the  MRDLL’s  finite  tracking  error;  this 
was  because  the  NCDLL  is  designed  assuming  no  multipath,  while  the  MRDLL  must  estimate  the 
received  signal  parameters  to  determine  the  strength  of  the  reflected  signal  or  if  the  reflected  signal 
is  present  at  all. 

The  final  simulation  investigated  MRDLL  performance  in  AWGN.  Results  indicated  that 
tracking  performance  predicted  by  linear  theory  was  applicable  only  when  the  received  SNR  ex¬ 
ceeded  a  given  threshold;  this  result  is  common  to  all  tracking  loop  analysis.  What  made  the 
MRDLL  SNR  results  unique  was  that  fact  that,  because  the  linear  tracking  region  depends  on  the 
relative  multipath  delay,  so  does  the  SNR  threshold;  as  one  might  expect,  the  greater  the  linear 
operating  region,  the  lower  the  required  SNR  threshold.  Results  demonstrated  that  this  type  of 
simulation  would  be  very  useful  in  determining  the  MRDLL’s  operating  specifications  for  different 
signal  parameters  (should  actual  GPS  parameters  be  used  in  the  simulation), 

5.4  Summary  and  Recommendations 

This  thesis  has  shown  that  The  MRDLL  has  tremendous  potential  for  GPS  application  where 
multipath  mitigation  is  a  concern.  This  was  demonstrated  by  both  simulation  and  analysis  in  which 
the  MRDLL  exhibited  significant  improvement  over  the  NCDLL  in  terms  of  steady-state  code-phase 
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tracking  error  with  a  single  reflection  present.  The  following  is  a  list  of  suggested  follow-on  work 
for  MRDLL  research  (and  research  on  GPS  multipath  mitigation  in  general): 

1.  Develop  a  MATLAB  (or  another  alternative)  simulation  model  of  the  MRDLL.  This  would 
make  simulations  more  efficient  than  they  were  in  SPW,  and  perhaps  allow  the  simulations 
to  run  at  GPS  parameters. 

2.  Design  a  detection/estimation  scheme  that  detects  the  presence  of  multipath,  and  allows  the 
GPS  receiver  to  switch  between  the  NCDLL  and  the  MRDLL  depending  on  the  received 
signal  environment.  This  could  optimize  code  phase  tracking  performance  since  the  NCDLL 
is  still  the  optimum  choice  for  GPS  code  phase  tracking  when  multipath  is  not  present. 

3.  Change  the  MCEU  estimation  scheme  to  allow  3  to  be  any  element  of  the  interval  (0,  1.5]. 
To  be  useful,  the  MCTL  would  also  have  to  be  modified  to  allow  mitigation  of  any  of  a 
continuum  of  possible  delays. 

4.  Generalize  the  mitigation  task  for  a  multiple-reflection  scenario. 

5.  Develop  a  carrier  loop  simulation  model  to  run  in  tandem  with  the  MRDLL  model. 

6.  Investigate  the  effects  of  code  and  frequency  Doppler  shift  on  the  MRDLL  in  the  presence  of 
multipath;  compare  both  ALC  and  non-ALC  loop  designs. 

7.  Investigate  MRDLL  performance  when  random  data  modulation  is  included  on  the  received 
GPS  signal.  The  data  can  then  be  demodulated  and  a  bit  error  rate  analysis  conducted. 

8.  Further  investigate  the  minimum  mean-square  error  estimation  scheme  suggested  by  Weill 
(14). 

9.  Consider  narrower  correlator  spacings  (A  <  1)  in  the  MCTL  (MRDLL). 
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Appendix  A,  Computer  Simulation  Models 

A.l  Received  GPS  signal  -  complex  envelope  representation 

Each  of  the  simulation  models  presented  in  this  appendix  are  based  on  the  complex  low- 
pass  envelope  of  the  received  GPS  multipath  signal.  The  complex  low-pass  approach  is  a  useful 
simulation  tool  because  it  decreases  the  required  sampling  frequency  and  simulation  time;  the 
approach  is  valid  as  long  as  the  carrier  frequency  is  much  greater  than  the  bandwidth  of  the 
bandpass  signal  (as  is  the  case  for  GPS)  (3).  To  determine  the  complex  envelope  representation, 
the  GPS  multipath  signal  (neglecting  data  modulation  and  noise)  is  written  as 

r{i)  —  y/2P  {aQc[t  -  ro)cos(a;ct  +  ^o)  +  aic{t  -  tq  -  aTc)  cos{wct  -f  ^i)}  (101) 

=  [ooc(t  -  +  aic{t  -  tq  -  ]■ 

=  Re  {r(t)e^‘*^^*} 

where  r{t)  is  the  complex  envelope  of  the  signal  given  by 

r{t)  =  VW  [aQc{t  ~  aic{t  -  To  -  aTc)e^^^]  ,  (102) 

Therefore,  as  shown  in  Figure  83,  the  signal  can  be  modeled  as  a  DS/SS  code  passing  through  a 
tapped  delay  line  with  complex  tap  weights  (7). 

For  the  MRDLL,  the  received  signal  must  be  converted  to  baseband.  Recall  from  Chapter 
III  that  this  conversion  is  accomplished  by  mixing  the  received  signal  (neglecting  data  and  noise) 
with  a  local  carrier  to  give 

rbb{t)  =  XQc{t  -  To) Xic{t  -  To  -  aTc)  (103) 

=  \/W{ao  cos{9o  -  63)0  {t  -  To)  -|-  ai  cos(^i  -  03)0  {t  -  tq  -  aTc)} 
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Figure  84.  SPW  GPS  Multipath  Signal  block  detail. 

where  Oz  is  an  arbitrary  phase  error.  Therefore,  the  baseband  conversion  is  modeled  by  multiplying 
r{t)  in  Equation  102  by  and  taking  the  real  part  to  give  Equation  103. 


A. 2  GPS  Multipath  Signal  model 

The  GPS  MULTIPATH  SIGNAL  block  detail  is  shown  in  Figure  84.  This  block  generates  the 
received  GPS  signal  according  to  Equation  102  and  the  channel  model  in  Figure  83;  the  direct-path 
and  multipath  components  of  the  signal  are  monitored  separately  via  the  ‘dp-Out’  and  ‘mp.out 
ports,  respectively.  The  GPS  MULTIPATH  SIGNAL  block  consists  of  two  main  components:  the 
PN  CODE  GENERATOR  and  the  MULTIPATH  CHANNEL.  Two  COMPLEX  DELAY  blocks 
(marked  'Z’  in  the  diagram)  are  also  included;  the  purpose  of  these  blocks  will  be  explained  in 
Section  A. 3.1. 

The  PN  CODE  GENERATOR  block  detail,  shown  in  Figure  85,  generates  a  DS/SS  code 
at  a  specified  chip  rate  (1/Tc)  and  code  period,  N.  The  desired  chip  rate  is  achieved  by  setting 
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Figure  85.  SPW  PN  Code  Generator  block  detail. 

the  IMPULSE  TRAIN  frequency  equal  to  1/Tc;  the  impulse  train  acts  as  a  control  input  to  the 
GENERAL  PN  SEQ  SOURCE  which  generates  a  majcimal-length  PN  code  of  Is  and  Os.  The  PN 
code  period  is  determined  by  setting  the  block’s  shift  register  order  according  to 

AT  =  2”*  -  1  (104) 

where  m  is  the  shift  register  order.  The  BPSK  MOD  block  then  converts  the  PN  code  to  a  binary 
waveform  of  Is  and  -Is  to  give  the  DS/SS  code,  c{t). 

The  MULTIPATH  CHANNEL  block,  shown  in  Figure  86,  implements  the  tapped  delay  line 
model  shown  in  Figure  83.  The  COMPLEX  VARIABLE  DELAY  blocks  allow  the  designer  to  vary 
the  values  of  tq  and  a,  while  the  G/P  MULT  blocks,  representing  the  complex  tap  weights,  allow 
the  designer  to  vary  the  attenuation  coefficients,  ai,  and  the  carrier  phases,  0^. 

A. 3  NCDLL  model 

The  NCDLL  block  detail,  shown  in  Figure  87,  is  based  on  the  NCDLL  block  diagram  and 
analysis  presented  in  Chapter  II.  Upon  receiving  the  multipath  signal,  the  NCDLL  block  returns 
two  outputs!  the  on-time  code  and  the  loop  filter  signal.  The  on-time  code  (  on-time  in  Figure  87) 
is  a  replica  of  the  incoming  direct-path  DS/SS  code  and  has  a  code  phase  equal  to  the  NCDLL’s 
estimate  of  the  direct-path  propagation  delay.  The  loop  filter  signal  (‘loop-flt’  in  Figure  87)  is  the 
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Figure  86.  SPW  Multipath  Channel  block  detail. 


lowpass  output  of  the  NCDLL  loop  filter  which  drives  the  VCC;  this  signal  is  monitored  to  ensure 
the  loop  remains  stable  and  to  determine  when  the  loop  reaches  steady-state. 

A,3J  Early/late  correlator  model.  Now,  we  will  examine  more  closely  the  individual 
components  of  the  NCDLL  block,  starting  with  the  early /late  correlators.  The  correlators  are 
modeled  as  multipliers  which  mix  the  incoming  GPS  signal  with  shifted  early  and  late  versions  of  the 
on-time  replica  code  generated  by  the  CODE  GENERATOR  block.  The  delay  block  preceding  the 
‘late’  correlator  ensures  that  the  signal  into  the  multiplier  lags  the  signal  into  the  ‘early’  multiplier 
by  one  chip  period,  Tc]  this  corresponds  to  an  early/late  offset  of  A  =  1. 

The  early/late  correlator  implementation  just  described  dictates  that  there  be  one  half-chip 
(T’c/2)  delay  between  the  CODE  GENERATOR  block  and  the  ‘on-time’  output;  this  delay  is 
accomplished  by  placing  a  delay  block  just  before  the  NCDLL  block’s  ‘on-time'  output.  The 
additional  unit  delay  block  immediately  after  the  CODE  GENERATOR  is  necessary  because  SPW 
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requires  at  least  one  sample  of  delay  in  feedback  loops.  The  presence  of  these  two  delay  blocks  is 
accounted  for  by  the  two  delay  blocks  in  the  GPS  MULTIPATH  SIGNAL  block  (Figure  84). 

A.3.2  Bandpass  and  lowpass  filter  models.  The  NCDLL’s  BPFs  are  modeled  as  LPFs 
since  we  are  using  complex  lowpass  signals.  Therefore,  these  filters  (along  with  the  LPFs  following 
the  squaring  operation)  are  modeled  using  SPW’s  COMPLEX  BUTTERWORTH  LOWPASS  HR 
FILTER  block;  this  block  allows  the  designer  to  specify  the  filter  passband  edge  (in  Hz)  and  the 
amount  of  attenuation  (in  dB)  at  the  passband  edge. 

A. 3.3  Squaring  circuit  model.  The  NCDLL’s  squaring  operation  is  modeled  by  the 
SQUARING  CIRCUIT  block.  Because  we  are  using  complex  lowpass  signals,  this  block  cannot 
simply  return  the  square  of  the  input  signal.  Referring  back  to  the  NCDLL  multipath  analysis  in 
Chapter  II,  the  NCDLL  early /late  correlator  output,  j/(f),  into  the  squaring  circuit  is  given  by 

y(t)  =  \/2p|aoi?c  0  cos(a;ot +^o)  +  aii?c  +  «±  0  cos(a;et  +  0i)|  ,  (105) 
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and  the  lowpass  component  of  the  squaring  circuit  output  (into  the  adder  circuit)  is 


j,ni)Lp  =  P  ±  i) +«!«*(<  + a  ±i)  (106) 

+  20^01  C03(^o  ±  2)  2)  }  ' 

The  complex  lowpass  envelope  of  y{t)  into  the  squaring  circuit  is 

y{t)  =  -v/^P  |aoi?c  ±  0  +  aiPc  ^<5  +  «  ±  ^)  |  •  (107) 

Multiplying  by  the  complex  conjugate,  y  *{t)  and  dividing  by  two,  we  see  that 

\yitW{t)  =  y\t)\^p  ■  (108) 

Therefore,  in  order  to  accurately  model  the  NCDLL  squaring  operation,  the  SQUARING  CIRCUIT 
block  multiplies  the  incoming  complex  lowpass  signal  by  its  complex  conjugate  and  divides  the  result 
by  two. 

AJ4  Loop  filter  model  The  loop  filter  is  modeled  by  SPW’s  FILTER  block  which  defines 
a  filter  according  to  its  difference  equation  coefficients.  As  seen  in  (3),  the  digital  realization  of  the 
active  lead*lag  filter  described  in  Equation  24  corresponds  to  the  difference  equation 

y  [n]  =  Aqx  [n]  +  Aix  [n  -  1]  -  Biy[n  -  1]  (109) 

where  x  [n]  and  y  [n]  are  the  discrete- time  representations  of  the  filter  input  and  output,  respectively. 
The  difference  equation  coefficients  are  given  by  (3) 
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Figure  88.  SPW  VCC  block  detail. 
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Figure  89.  Illustration  of  VCC  operation  in  SPW. 
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where  is  the  sampling  period  in  seconds;  the  FILTER  block  allows  the  designer  to  specify  these 
coefficients  according  for  the  desired  values  of  ri  and  T2. 

A. 3. 5  VCC  model.  The  VCC  block  adjusts  the  frequency  of  the  on-time  code  so  that 
its  phase  matches  the  NCDLL’s  estimate  of  direct-path  phase.  The  VCC  block  detail  is  shown  in 
Figure  88  and  its  operation  is  illustrated  in  Figure  89;  this  implementation  is  based  on  a  design 
suggested  in  (3).  The  VCC  output  serves  the  same  purpose  in  the  CODE  GENERATOR  block 
as  the  IMPULSE  TRAIN  does  in  the  PN  CODE  block  (see  Figure  85);  the  difference  is  that  the 
VCC’s  impulse  train  output  changes  in  frequency  according  to  the  VCO  output  as  shown  in  Figure 
89. 
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The  VCC  block  allows  the  designer  to  specify  the  VCO  quiescent  frequency  in  Hz  and  the 
VCO  constant  in  Hz/volt.  The  quiescent  frequency  of  the  VCO  (i.e.,  the  output  frequency  with  a 
zero-volt  input)  is  set  equal  to  the  chip  rate,  l/Tc  Hz.  Therefore  (with  no  noise  and/or  multipath), 
if  the  on-time  replica  code  is  in  phase  with  the  received  signal  code,  the  early/late  correlator  outputs 
will  be  equal.  This  will  result  in  a  loop  filter  output  equal  to  zero  and  will  cause  the  replica  code  to 
remain  at  1/Tc  Hz.  However,  if  the  received  code  lags  the  on-time  code,  the  received  code  will  more 
closely  correlate  with  the  late  replica;  therefore,  the  loop  filter  output  will  be  positive.  Similarly, 
if  the  received  code  leads  the  on-time  replica,  the  loop  filter  output  will  be  negative.  Therefore, 
the  VCO  constant  is  specified  in  such  a  way  that  the  VCC  frequency  increases  with  a  negative 
input  voltage  and  decreases  with  a  positive  input  voltage.  The  actual  value  of  the  VCO  constant 
is  chosen  according  to  the  desired  dynamic  performance  of  the  loop. 


A. 4  MCTL  model 

The  MCTL  block  detail,  shown  in  Figure  90,  is  based  on  the  MRDLL  block  diagram  and 
MCTL  analysis  presented  in  Chapter  III.  Like  the  NCDLL  block,  the  MCTL  block  generates  an 
on-time  replica  of  the  direct-path  code  and  a  loop  filter  signal  for  monitoring  purposes.  Unlike 
the  NCDLL  block,  however,  the  MCTL  block  receives  signal  parameter  estimates  from  the  MCEU 
block  in  addition  to  the  received  baseband  signal. 

The  MCTL  model  has  two  main  arms,  each  consisting  of  an  EARLY-LATE  CORRELATOR 
block  followed  by  a  gain/phase  multiplier.  The  top  arm  represents  the  MCTL’s  direct-path  arm, 
and  the  lower  arm  represents  the  enabled  multipath  arm  corresponding  to  the  MCEU’s  estimate  of 
the  multipath  delay  coefficient,  S.  Both  arms’  EARLY-LATE  CORRELATOR  blocks  (see  Figure 
91)  operate  in  the  same  manner  as  the  NCDLL’s  early /late  correlator  arms  (except  for  the  squaring 
operation).  In  the  multipath  arm,  a  delay  block  delays  the  MCTL’s  on-time  replica  by  aTc  seconds 
prior  to  early/late  correlation.  Therefore,  each  arm’s  early/late  correlator  output  is  the  correlation 
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Figure  90,  SPW  MCTL  block  detail, 
between  the  received  signal  and  a  shifted  replica  code  given  by 

ca  (<)|<jp  =  c(i  -  To  -  Tc/2)  -  c(i  -  To +Tc/2)  (111) 

for  the  direct-path  arm  and 

{t)\mp  =  c{t-To  -STc-  Tc/2)  -c{t-To-  STc  +  1^/2)  (112) 

for  the  enabled  multipath  arm.  The  early/late  correlator  outputs  then  enter  the  gain/phase  mul¬ 
tipliers,  where  they  are  mixed  with  the  appropriate  MCEU  signal  estimates.  If  the  MCEU  output 
is  S  =  0  (i,e.,  no  multipath  present),  a  switch  disables  the  lower  arm. 

The  ALC  block  operates  in  conjunction  with  the  SPW  FILTER  block  to  comprise  the  overall 
ALC  model  (refer  to  Figure  19  for  the  ALC  block  diagram).  The  ALC  block  receives  signal  estimates 
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Figure  91.  SPW  Early-Late  Correlator  block  detail. 


from  the  MCEU  block  and  adjusts  the  loop  filter  gain,  K^,  as  described  in  Chapter  III;  the  adaptive 

gam  threshold,  is  also  set  within  this  block.  In  the  FILTER  block,  the  difference  equation 

coefficients  in  Equation  110  are  chosen  so  that  the  loop  filter  transfer  function  is  described  in 
Equation  61 

A. 5  MCEU  model  (MAT LAB) 

A. 5.1  Description.  As  seen  in  Chapter  III,  the  MCEU  relies  on  matrix  operations  to 
compute  the  multipath  signal  estimates  sent  to  the  MCTL.  The  MATLAB  software  package  is 
optimized  for  matrix  operations;  therefore,  a  MATLAB  model  of  the  MCEU  was  developed  to 
perform  preliminary  simulation  and  analysis  prior  to  SPW  implementation. 

The  MCEU  MATLAB  simulation  model  is  based  on  the  MCEU  block  diagram  in  Figure 
25  and  the  analysis  presented  in  Chapter  III.  The  MATLAB  code  for  this  simulation  model  is 
contained  m  the  mcea.m  m-file  m  Section  A.5.2.  Running  mceu.m  represents  one  sample  of  the 
MCEU’s  M  signal  parameter  estimates  (x,,  and  3^)  and  their  corresponding  error  estimates  E,. 
The  model  assumes  a  received  baseband  signal  as  described  by  Equation  103,  zero  tracking  error 
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(i,e.,  To  =  7^),  and  large  code  period,  iV,  so  that  the  code  autocorrelation  approximation  in  Equation 
5  holds.  Also,  we  choose  the  number  of  correlator  arms  as  M  =  16  with  replica  codes  spaced  O.lTc 
seconds  apart  (i.e.,  Pk  =  O.lfc). 

Care  must  be  taken  when  modeling  the  correlator  noise  outputs  in  MATLAB.  Recall  that 
the  correlator  output  noise  samples  are  not  uncorrelated;  instead,  they  are  represented  as  a  16  x  1 
random  vector,  v,  distributed  as  N  :  [O,  o-^Cv],  where  al  is  the  variance  of  v  and  Cv  is  the  noise 
covariance  matrix  defined  in  Equation  75.  MATLAB  can  generate  a  similar  random  vector,  w,  but 
with  distribution  N  :  [0,1],  where  I  is  the  identity  matrix.  Therefore,  it  is  necessary  to  transform 
the  MATLAB-generated  white  noise  vector,  w,  to  the  colored  noise  vector,  v.  The  whit2col.m 
m-file  (Section  A.5.3)  accomplishes  this  transformation  according  to  (6) 

V  =  ctvL’^w  (113) 

where  cr^  is  the  noise  standard  deviation  and  L  is  the  Cholesky  factorization  of  Cv 
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A,b.2  The  mceu  M-file. 


y.MCEU 

•/.Multiple-Correlator  Estimation  Unit  (MCEU)  estimator  output 

y. 

y,[E,  X,  ALPHA_ML,  X_ML]  =  MCEU(X0,  XI,  ALPHA,  CNo,  B)  returns  the  MCEU  estimator 
output 

Xfor  received  signal  parameters,  XO,  XI,  and  multipath  delay  coefficient,  alpha, 
yihe  remaining  two  input  aorguments  are  the  received  direct-path  carrier-to-noise 
density 

•/ratio  (CNo  in  dB-Hz)  calculated  for  received  carrier  power  P=l/2,  and  the  MCEU 
•/integrator  one-sided  lowpass  bandwidth  (B  in  Hz)  (setting  B=0  gives  no-noise 
response) . 

•/ 

•/The  output  arguments  are  defined  as  follows: 

y, 

•/E  =  a  vector  whos  elements  are  the  16  candidate  error  measurements,  E_k;  (the 
MCEU 

•/chooses  as  its  ML  estimate  the  candidate  signal  parameter  estimate 
corresponding  to 
•/the  minimum  E^k)  . 

% 

•/X  =  a  2  X  16  matrix  whos  columns  correspond  to  the  16  candidate  signal 
parameter 

Xand  rows  correspond  to  xO  (1st  row)  and  xl  (2nd  row) 

1 
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7,ALPHA_ML  =  The  MCEU’s  ML  alpha  estimate 
•/o 

'/.X^ML  =  The  MCEU^s  ML  signal  parameter  estimate  vector 

t 

VoThis  fxinction  uses  the  WHIT2C0L  m-file  to  generate  the  correlator  output  noise 
•/.vector,  and  the  CODECORR  m-file  to  compute  the  code  autocorrelation  function. 

I 

•/.Written  by  Mark  C.  Laxton,  1  May  96.  Last  updated:  15  Nov  96. 

function  [E , X ,  alpha_ml  ,x_ml]  =mceu(xO  ,xl ,  alpha, CNo  ,B) 

beta=0:0.1:1.5;  '/.Set  up  delay-spacing  coefficients,  beta 
alpha_k=beta;  ‘/.Set  alpha  estimates,  alpha_k,  equal  to  beta 

v=whit2col(beta,CNo,B) ;  •/.Generate  correlator  noise  vector,  v 
v=v' ; 

•/.Define  correlation  matrix,  C 
w=codecorr(beta-beta(l)) ; 

C=toeplitz(w) ; 

•/.Form  correlator  measurement  vector,  R 
R=codecorr (beta) ♦xO+codecorr (alpha-beta) *xl+v; 

•/.Form  candidate  estimates  for  xO  cind  xl  according  to  MCEU  estimation  algorithm 
for  k=0 : length (alpha_k)-l 
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if  k==0 


H=(codecorr(beta)) ^ ; 

x0_0=inv (H  ^  * inv (C) *H) *H ’ *inv (C ) *R ’ ; 

xl_0=0; 

X(:,k+l)  =  [x0_0  xl^O]^ 
else 

H( : , l)=(codecorr (beta) ) ^ ; 

H( : ,2)=(codecorr(alpha_k(k+l)-beta)) ^ ; 

X(: ,k+l)=inv(H'*inv(C)*H)*H'*inv(C)*R' ; 
end 
end 

7oClioose  ML  estimate  by  computing  and  choosing  estimates  corresponding  to  minimum 
E. 

for  k=0 :length(beta)-l 
H( : , l)=(codecorr(beta) ) ^ ; 

H( : ,2)=(codecorr(alpha_k(k+l)-beta)) ’ ; 

R.est( : ,k+l)=H*X( : ,k+l) ; 
end 

R_chkmx=R^  ♦ones  ( 1 ,  length  (bet  a)  )  ; 

e=R_est-R_chkmx ; 

for  k=0 : length (bet a) -1 

E(k+l)=e( : ,k+l) ^ *inv(C)*e( : ,k+l) ; 

end 


m=find(E==min(E)) ; 
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A. 5. 3  The  whit2col  M-file. 


y,WHIT2C0L 

•/.Transform  white  noise  to  colored  noise. 

•/• 

•/.(This  function  was  written  to  be  used  in  conjunction  with 
•/.the  MCEU  m-f  ile)  . 

7. 

*/.Z  =  WHIT2C0L(BETA,CNo,B)  transforms  a  white  noise  vector  distributed  as 
7.N[0,var*I]  to  a  noise  vector,  Z,  distributed  as  N[0,var*R]  where  R  is  the 
•/.correlator  matrix  corresponding  to  the  noise  outputs  of  the  MCEU  correlators 
•/.with  delay  spacing,  BETA.  The  integrator  bandwidth  of  the  BPFs  is  B  (in  Hz), 
•/.and  the  received  carrier-to-noise  density  ration  is  CNo  dB“Hz  (assumes  signal 
•/.power,  P=l/2). 

7. 

•/.Written  by  Mark  C.  Laxton,  1  May  1996.  Last  updated:  15  Nov  96 

function  z=whit2col (beta, CNo ,B) 

•/.Create  white  noise  vector,  w 
k=length(beta) ; 
w=randn(k,l) ; 

•/.define  correlation  matrix,  R 
v=codecorr (beta~beta(l) ) ; 

R=toeplitz(v) ; 
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y.Use  Cholesky  factorization  to  determine  transformation  matrix  A 


A=chol(R); 


•/transform  w  to  output  noise  vector,  z 
z=A^*w; 


•/multiply  by  sqrt  of  varicince,  sig2 
No=l/(2*10*(CNo/10)); 
sig2=2*B*No ; 
z=sqrt(sig2)*z; 
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A. 5. 4  The  codecoTT  M-file. 


7.C0DEC0RR 

•/.Compute  DS/SS  code  autocorrelation  function  for  large  code  period 

7. 

•/,Rc  =  CODECORR(x)  returns  the  autocorrelation  function,  Rc,  of  a  DS/SS 
•/^code  having  a  code  period,  N  »  1.  The  value  of  Rc  is  approximated  by 
•/,Rc(x)  =  1-1x1  for  |xl<=l  and  Rc(x)=0  elsewhere. 

•/. 

•/.Written  by  Mark  C.  Laxton,  1  May  1996.  Last  updated;  15  Nov  96 

function  Rc=codecorr(x) 

for  k=l :length(x) 

if  abs(x(k))<=l 

Rc(k)=l-abs(x(k)) : 

else 

Rc(k)=0; 

end 

end 
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Figure  92.  SPW  MCEU  model. 


A.6  MCEU  model  (SPW) 

The  SPW  MCEU  model  is  shown  in  Figure  92  and  consists  of  two  main  components:  the 
multiple  correlator  bank  (CORRBANK)  block  and  the  maximum-likelihood  estimator  (MCEU -ML) 
block.  The  CORRBANK  block  generates  M  =  16  correlator  measurements  that  are  fed  into  the 
MCEU-ML  block;  each  correlator  is  modeled  using  the  variable-delay  correlator  (VDC)  block  shown 
in  Figure  93.  Each  VDC  correlates  the  received  baseband  signal  with  a  delayed  version  of  the  on- 
time  replica  code  from  the  MCTL  block;  the  fcth  VDC  has  a  code  replica  delay  spacing  coefficient, 
Pk  =  O.lfc.  The  MCEU-ML  block  receives  the  CORRBANK  correlator  measurements,  computes 
the  M  =  16  signal  parameter  estimates,  and  sends  the  ML  estimate  to  the  MCTL  block  as  described 
in  Chapter  III. 
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code 


Figure  93.  SPW  Variable-Delay  Correlator  (VDC)  block  detail. 


Appendix  B.  Data  Processing  M- files 


B.l  The  mrdlest  M-file 

*/,HRDLEST  Compute  and  plot  means  and  variances  of  MCEU  ML  output  estimates 

•/. 

'/.EST  =  MRDLEST(START,STOP,XOFNM,XlFNM,ALPHAFNM,S_FREq)  computes  the  sample  means 
Viand  Vctriances  of  three  vectors  (‘xO’,  ‘xl’,  and  ‘alpha’)  which  contain  samples 
•/.of  the  corresponding  MCEU  ML  estimator  outputs;  the  sample  frequency  is  defined 
*/,by  the  S_FREQ  input  argument.  The  three  vectors  are  contained  in  mat-files  named, 
•//XOFNM’,  ‘XIFNM’,  and  ‘ALPHAFNM’.  Each  mat-file  is  assumed  to  have  come  from 
'/.an  SPW  MATLAB  SINK  which  has  assigned  the  output  data  to  a  vector,  ‘Y’ ,  and 
•/.has  also  returned  an  erronneous  ‘sampfreq’  value  (based  on  SPW  v.3.5,  1996).  This 
•/, m-file  renames  each  Y  vector  to  xO,  xl,  or  alpha  as  appropriate,  and  deletes  the 
^eiToneous  sampfreq  variable .  The  range  of  data  samples  used  to  compute  the  sample 
'/mieans  and  variances  are  defined  by  the  START  and  STOP  input  arguments;  these  values 
•/.specify  the  first  and  last  indices  of  the  xO,  xl,  and  alpha  vectors  used  in  the 
'/.mean  and  variance  computation.  The  output  matrix,  EST,  is  a  2  X  3  matrix  whose 
'/.columns  correspond  to  xO,  xl,  and  alpha  (in  that  order),  and  whose  1st  row  is  the 
'/.sample  mean  and  2nd  row  is  the  sample  variance. 

•/. 

•/.EST  =  MRDLEST (..  ..LABEL,  XFIG,  AFIG)  plots  the  MCEU  outputs  vs  time  and  sends  the 
•/.results  to  a  print  file  in  eps  format .  The  LABEL  argument  is  a  string  containing 
•/.the  title  of  the  plot;  XFIG  and  AFIG  are  strings  which  correspond  to  the  print 
'/.filenames  of  each  plot. 

•/. 

•/.Written  by  Mark  C.  Laxton,  1  May  96.  Last  updated:  15  Nov  96. 
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function  EST=iiirdlest(start , stop, xOfnm,x If nin,alphafniii,s _freq, label ,xfig,afig) 

•/.load  matfiles,  rename  Y  vectors,  clear  sampfreq  value 
evaKC’load  ^  xOfnm  *  -mat’]); 
clear  sampfreq 
xO=Y; 

eval([’load  ’  xlfnm  ’  -mat’]); 

clear  sampfreq 

xl=Y; 

eval([’load  ’  alphafnm  ’  -mat’]); 
clear  sampfreq 
alpha=Y ; 
clear  Y 

Voif  plot  options  are  selected,  plot  vs  time  and  send  to  printfile 

if  nargin>6 

t=0 : length(xO) -1 ; 

time=t/s_freq; 

f igure, plot (time, xO) 

hold  on 

plot (time, xl, ’ — ’) 
hold  off 

xlabeK’time  (sec)’) 
ylabeK’MCEU  output’) 
grid 
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s=iiiim2str(label) ; 


title(s) 

h=legend(^-' , ’xO  estimate \ 'xl  estimate \0); 
axes(h) 

eval(  [’print  -deps  -epsi’  xfig]); 

figure ,plot (time , alpha) 

grid 

xlabeK’time  (sec)’) 
ylabeK’MCEU  output’) 
title(s) 

h=legend( ’ alpha  estimate ’ ) ; 
axes(h) 

eval( [’print  -deps  -epsi’  af ig] ) ; 
end 

y, determine  mean  and  variance  of  estimates 
xOmn=me an (xO (start :stop)) ; 
xlmn=mean(xl (start : stop) ) ; 
alphainn=mean(  alpha  (start :  stop)  )  ; 
xOvcir=cov(xO( start :  stop) ) ; 
xlvcir=cov(xl  (start :  stop) ) ; 
alphavcLr=cov ( alpha (st cirt :  stop)  )  ; 

EST=[xOiTin  xlmn  alphainn ;  xOvar  xlvar  alphavar]  ; 
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B.2  The  leadedge  M-file 


•/.LEADEDGE 

'/^Normalized  rms  code  phase  tracking  error  estimate . 

% 

%  RMS_ERR  =  LEADEDGECDP,  OT,  CHP.RATE,  S^FREQ)  calculates  the  rms 
%  code  phase  timing  error  (in  chips)  between  a  direct-path  version  of  a 
%  BPSK -modulated  (l^s  and  -I’s)  direct-sequence  spread- spectrum  (DS/SS)  code 
•/,  and  and‘ on-time^  replica  generated  by  a  code  tracking  loop.  This  is  done  by 
•/,  comparing  the  leading  edges  between  the  two  codes.  The  direct-path  and 
'/o  on-time  codes  are  input  as  mat-file  filenames  (^DP^  and  'OTO  for  a 
%  specified  code  chip  rate  (CHP_RATE)  in  Hz  and  simulation  sampling  freq 
%  (^S^FREQO  in  Hz. 

7o 

7,  RMS_ERR  =  LEADEDGECDP,  OT,  CHP3ATE,  S^FREQ,  LOCK)  will  display  a 

•/,  "loss  of  lock"  message  if  the  code  phase  error  between  any  pair  of  leading 

'/,  edges  exceeds  the  threshold  (in  chips)  specified  by  LOCK. 

7« 

%  RMS^ERR  =  LEADEDGECDP,  OT,  CHP3ATE,  S_FREQ,  LOCK,  LIN)  will  display  an 
•/,  additional  "out  of  linear  region"  message  if  the  code  phase  error  between 
t  any  pair  of  leading  edges  exceeds  the  threshold  (in  chips)  specified  by  LIN. 
% 


I  CAUTION:  Care  must  be  taken  to  ensure  that  the  first  pair  of  leading  edges 
'/a  correspond  to  the  same  point  in  the  DS/SS  code  sequence. 


/#  /o  /#  /#  /o  h  h  h  h  h  h  h  h  /o  h  /o  h  h  h  h  h  h  /#  h  /o  h  h  h  h  h  h  h  h  h  h  /o  h  h  h  /o  /« h  h  h  h  h  h  h  h  h  h  h  h  h  /o  h  /o  h  /o  /o  /a  /•  /•  /o  /o 


'/o  Written  by  Mark  C.  Laxton,  1  May  96.  Last  updated:  15  Nov  96. 


function  rms_.err=leadedge(dp ,  ot ,  chp.rate ,  s„f req,  lock ,  lin) 


•/.Find  samples  where  each  code  equals  -1;  put  these  samples  into  two  vectors, 

•//fdp^  and  ^fot^ 
fdp=f ind(dp==-l) ; 
f ot=f ind(ot==-l) ; 

•/.Determine  which  samples  correspond  to  leading  edges  (-1  to  1  transitions) 

m=l; 

p=i; 

for  k=l :length(fdp)-l 

if  fdp(k)+l"=fdp(k+l)  •/.Check  if  each  sample  number  in  ’fdp’  is  followed  by  the  next 
•/.highest  integer.  If  not,  that  sample  is  a  leading  edge. 

ledp(m)=fdp(k) ;  '/.Form  vector  ’ledp'  of  sample  numbers  corresponding  to  the 

y.leading  edges  of  the  direct -path  code . 

m=m+l; 

end 

end 

for  k=l:length(fot)-l  •/.Repeat  for  on-time  code,  forming  ^leot'  vector. 

if  fot(k)+l"=f ot (k+1 ) 

leot(p)=f ot (k) ; 

p=p+l; 

end 

end 
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VoMake  the  number  of  leading  edges  for  both  codes  the  same . 

if  length(ledp)>length(leot) 

ledp=ledp(l:leiigth(ledp)-l)  ; 

elseif  length(ledp)<length(leot) 

leot=leot ( 1 : length(leot ) -1 ) ; 

end 

•/oDetermine  number  of  samples  between  each  pair  of  leading  edges,  place  these  values  into 
•/oa  vector,  ^sampdiff^. 
sampdif  f =leot-“ledp ; 

/^Determine  ’loss  of  lock’,  and  ’out  of  linear  region’  outputs, 
if  nargin>4 

loselock=find(abs(sampdiff  )>(lock*s_freq/chp^rate))  ;  ’/oDetermine  the  samples  where  the 
•/o’loss  of  lock’  threshold  is 
Vo  is  exceeded, 
if  loselock"'=[] 

disp([’  ’;’loss  of  lock’])  7Jf  the  threshold  has  been  exceeded  at  any  two  pair 

'/oOf  leading  edges,  display  ’loss  of  lock’. 

end 

end 

VoRepeat  similar  steps  for  ’out  of  linear  region’ 
if  nargin>5 

out_lin=find(abs (sampdif f )>(lin*s«freq/chp_rate) ) ; 
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if  out.lin"'=[3 


dispCC^  ^;’out  of  linear  region^]) 

end 

end 

y, Calculate  nhs  error  (in  chips) 

rms^err=chp_rat€*sqrt(mean(sampdiff .  "2))/s^freq; 
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Appendix  C.  Simulation  #4-  MCEU  ML  Estimator  Variance  in  AWGN 


MRDLL  SNR  simulation,  a  =0.4 


Figure  94.  Simulation  results  showing  the  variance  of  the  MCEU  ML  estimate,  xq  for  different 
values  of  effective  MCTL  loop  SNR,  peff  (in  dB)  when  a  =  0,4. 


MRDLL  SNR  simulation,  a  =0.4 


Figure  95.  Simulation  results  showing  the  variance  of  the  MCEU  ML  estimate,  Xi  for  different 
values  of  effective  MCTL  loop  SNR,  peff  (in  dB)  when  a  =  0.4. 
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MRDLL  SNR  simulation,  a  =0.5 


Figure  98. 


Figure  99, 


Simulation  results  showing  the  variance  of  the  MCEU  ML  estimate,  x\  for  different 
values  of  effective  MCTL  loop  SNR,  peff  (in  dB)  when  o;  =  0,5. 


MRDLL  SNR  simulation,  a  =0.5 


Simulation  results  showing  the  variance  of  the  MCEU  ML  estimate,  a  for  different 
values  of  effective  MCTL  loop  SNR,  peff  (in  dB)  when  a  =  0,5. 


MRDLL  SNR  simulation,  a  *0.7 


Figure  100.  Simulation  results  showing  the  variance  of  the  MCEU  ML  estimate,  Xq  for  different 
values  of  effective  MCTL  loop  SNR,  (in  dB)  when  a  =  0.7. 


MRDLL  SNR  simulation,  a  =0.7 


Figure  101.  Simulation  results  showing  the  variance  of  the  MCEU  ML  estimate,  xi  for  different 
values  of  effective  MCTL  loop  SNR,  peff  (in  dB)  when  a  —  0.7. 
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MRDLL  SNR  simulation,  a  =0.7 


Figure  102. 


Figure  103. 


Simulation  results  showing  the  variance  of  the  MCEU  ML  estimate, 
values  of  effective  MCTL  loop  SNR,  Peff  when  a  —  0.7. 


MRDLL  SNR  simulation,  a  =1.1 


Simulation  results  showing  the  variance  of  the  MCEU  ML  estimate, 
values  of  effective  MCTL  loop  SNR,  peff  (in  clB)  when  a  =  1.1. 


for  different 


for  different 


MRDLL  SNR  simulation,  a  =»1 .1 


Figure  104.  Simulation  results  showing  the  variance  of  the  MCEU  ML  estimate,  xi  for  different 
values  of  effective  MCTL  loop  SNR,  peff  (in  dB)  when  a  =  1.1. 


MRDLL  SNR  simulation,  a  =1 .1 


Figure  105,  Simulation  results  showing  the  variance  of  the  MCEU  ML  estimate,  a  for  different 
values  of  effective  MCTL  loop  SNR,  pe/f  (in  dB)  when  a  =  1.1. 


MRDLL  SNR  simulation,  a  =1 .5 


Figure  106.  Simulation  results  showing  the  variance  of  the  MCEU  ML  estimate,  xq  for  different 
values  of  effective  MCTL  loop  SNR,  peff  (in  dB)  when  a  =  1.5. 


MRDLL  SNR  simulation,  a  -1 .5 


Figure  107.  Simulation  results  showing  the  variance  of  the  MCEU  ML  estimate,  xi  for  different 
values  of  effective  MCTL  loop  SNR,  peff  (in  dB)  when  a  =  1.5. 
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MRDLL  SNR  simulation,  a  >1.5 


Figure  108.  Simulation  results  showing  the  variance  of  the  MCEU  ML  estimate,  a  for  different 
values  of  effective  MCTL  loop  SNR,  pe/f  (in  dB)  when  a  =  1.5. 
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